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Summary The theory of distillation is applied to approximate a multicomponent mixture by 
a simpler one. In both the top section and the bottom section, two substitute components can 
be found which provide the best approximation of the specified components by a binary mixture 
for a given moderately large number of trays (e.g. 10 or 20). The distillation can then be calcu- 
lated by the common methods applicable to a binary mixture. A third component is introduced 
if higher numbers of trays are considered but the calculation method still remains essentially 
that of a binary mixture. The unspecified components are taken into account by a simple 
correction method. Graphical methods are given to improve the latter correction and to calculate 
liquid or vapour rate of any actual component on an arbitrary tray. A constant ratio of 
volatilities is assumed. 


Résumé—L auteur applique la théorie de la distillation pour ramener le cas d'un mélange a 
plusieurs composants au cas d'un mélange binaire. Dans une section en téte et en queue de 
colonne, il peut trouver deux composants “ fictifs " qui permettent la meilleure approximation 
des composants * actuels * par un mélange binaire et cela pour un nombre de plateaux relative- 
ment grand (10 & 20). Les caleuls de la distillation s’effectuent alors par les procédés classiques 
applicables & un mélange binaire. Pour un plus grand nombre de plateaux, il introduit un 
troisiéme composant mais le calcul reste encore, essentiellement celui d'un mélange binaire. 
Par une méthode simple de correction, il tient compte des composants * négligés."" Par des 
méthodes graphiques, il peut ameéliorer cette correction et calculer le taux liquide et vapeur de 
chaque composant réel sur un plateau arbitraire. Il admet un rapport constant des volatilités. 


INTRODUCTION components replace the specified components of 

Several methods have been suggested for the the actual distillation. |The unspecified compo- 
approximation of a multicomponent mixture to a nents are separately taken into account [3, 4, 5}. 
binary one with a view to simplify the calculations The method is based upon an exact theory of 
In these methods, however, a certain arbitrariness. multicomponent distillation developed by one of 
is admittedly involved in the choice of the two the authors (6). In this present article a refined 
substitute components and it is diflicult to esti- method will be discussed resting upon the same 
mate the accuracy of the approximation [1, 2]. principles but for a great number of trays. Three 
The present authors have given a method in components instead of two are used to approxi- 
which the absorption factors on the trays in the mate the specified components in the actual 
column for the actual multicomponent distillation mixture, but the calculation method remains 
are approximated by the absorption factors for that for a binary mixture. A constant relative 
another distillation in which two substitute volatility in each section is assumed and a 
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constant reflux will be used in this article to 
simplify the notation. The latter restriction is 
however not essential since it has been shown that 
a distillation with variable reflux can be reduced 


to a distillation in which the reflux is constant [7]. 


TueoreticaL Basis or THe Meruop 


In a previous communication an exact equation 
has been derived from which the absorption 
factors of an arbitrarily chosen reference com- 
ponent can be calculated on each tray in the 
column. The absorption factor of a component 
ion a tray with index m is 


A L_/V.K, (1) 


Here A, ,, is the absorption factor, L, and V,, 
the total molar liquid and vapour flow from the 
mth tray respectively and A, ,, is the equilibrium 
constant for component ¢ at the mth tray.* If 
K,, is written for the equilibrium constant of the 
has the 


reference one 


assumption of constant relative volatility : 


K a,K, so A 


component owing to 


The stripping factor is : 


In these formulae L, and V,, are constant 
and L,, V, is written henceforth in the stripping 
section and L,, V, in the top-section. 

The product the 
stripping section of the reference component, 


Ap (m), is now introduced, It is by definition : 


A, Ay... 


t 


of absorption factors in 


Ap(m) = A A,, and Ap (0) = 1, (3) 


a 
if the plates are numbered upwards starting 
with m = 0 for the reboiler. 

The theory then leads to the following summa- 
tion equation from which Ap(m) can be calculated 


[1]: , 


L, Ap(m) = 2; = 


Z BaP Aplk) +2 


k=0 


' Vs ot 


ApM+1) V, V7" 


‘ m-(M +1) (4) 


Vw, 1 x, 


The notation 5.’ means a summation over all 


: 
components which have a specified residue rate 


* A list of notations is given at the end of the paper. 
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B, and &," is a summation over the volatile 
components the residue rate of which must be 
very small but is unspecified, V, y., is_ the 
vapour flow of such a component from the feed 
tray. The index M + 1 is used for the feed tray 
if it is considered as a tray of the stripping 
1 is used if it is considered 
The solution of 


section whereas N 
as a tray of the top section. 
equation (4) is 


Ap(m) Ap,(m) Ap(M 


tees 1) (5) 
The solution of equation (4) if all unspecified 


components are omitted is denoted here by 


Ap,(m) in this formula and C, is : 


I, Vu N41 
y 
t 
L, 


x," u (m—k) 
ke@ 
in which the function u (m k) is an abbrevia- 
tion for 2, B,a* Thus the solution of the 
equation in which the unspecified components 
are neglected must first be found. 


Specified Components 
Neglecting the unspecified components in equation 
(4) and assuming that the reboil vapour has the 
same composition as the residue or V; , = Ry, B; 
one obtains :* 


kom 
L, Ap, (m) z 
k 


u(m 
0 


k) Ap, (k) + Rygu(m) (7) 


the has 


equation 


The individuality of components 
disappeared from (7). The 
absorption factor is obtained with different mix- 
tures provided that the function u(m) remains 
the The specified components of the 
multicomponent therefore, be 
approximated by a binary mixture in such a way 
that the function : 


same 


same. 


system must, 


u*(m) = By, a," + By, a,” (8) 


is the best approximation of u(m). The substitute 


* A reboil vapour in equilibrium with the residue leads 
to Vig = So Bpa; and then 2, V; 4 a;" = Sg u(m+1) 
if S, is the stripping factor of the reference component in 
the reboiler. Thus equation (7) remains of the same type 
and its solution is obtained in the same way as above. 
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components have a relative volatility z,, and x,, 


b 
respectively, with respect to the veheunee com- 
ponent in the multicomponent mixture. Their 
residue rates are B,, and B,,. 

It is not necessary that u*(m) approximates 
In fact if the term u(0) Ap,(m) 
occurring in the summation in equation (7) is 
transferred to the left-hand further 
u*(O)Ap,(m) — (B,, B,,) Ap, (m) is added at 
both sides and finally u*(m) is substituted instead 


u(m) for m 0. 


side. if 


of u(m) one obtains an equation that approximates 


equation (7) from m l on: 

{L, — u(O) + u*(O)} Api) 
kom 
» u*(m kr) Ap, (k) R,, u*(m) (9) 
kh 0 


Let us now compare this equation with that 
for the binary mixture consisting of the com- 
ponents 1b and 2b with the residue rates B,, and 
and a liquid molar flow of 


u*(O). 


By, respectively 
L,* = L, + u(0) 
equation is according to equation (7) : 


Its summation 


hom 
L,* Ap,(m)= 2 u*(m—k)Ap,(k) + Ry*u*(m) (10) 
k 0 
In this equation R,* V,* u*(O) 
{L,* u*(0)} u*(O). It is identical to equation 


(9) except for the change of R, into R,*. It 
can be easily shown that the solution of equation 
(10) is equal to k,,* R, times the solution of 
equation (9) (vide. equation (8), p. 940 in [1}). 
Therefore the absorption factor of the reference 
component on the mth tray which is A,,, 
Ap, (m) Ap, (m 1) from 
equation (10) if m = 2 or higher since the factor 
R,* R, is divided out. We 
A, Ap, (1) Ap, (0) and since Ap, (0) 
both equation (9) and equation (10) we find an 
A, from equation (10) which is R,* R, times 
that of equation (9). 

Summarising we have shown that the product 


can be calculated 


have. however. 
1 for 


of absorption factors for the reference com- 
ponent calculated from the binary mixture to 
which equation (10) applies is equal to (R,* R,) 
Ap,(m) in which Ap,(m) is the solution of 
equation (9) which is an approximation of the 
absorption factor product of the multicomponent 
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mixture omitting the unspecified components 
(equation (7)). 

The degree of accuracy with which the solution 
of equation (9) approximates the absorption 
factor product of the multicomponent mixture 
depends on how accurately u* (m) approximates 
u(m) on the average since negative and positive 
affect the 
By applying the method of successive 


deviations solution in an opposite 
direction. 
approximations any desired degree of accuracy 
can be obtained. 

In the previous article the substitute com- 
ponents have been so chosen as to approximate 


u (m) for all values of m from m 1 on upwards 


[4]. 


u (m) becomes preponderant for large values of m 


Since the contribution of the light key to 


one of the substitute components had to be the 
light key itself thus leaving only one component 
which can be adapted for the approximation of 
u(m) by u* (im). 

If. however, the substitute components are so 
chosen that u* (m) approximates u(m) only for 
a finite range of values both components can be 
adapted. In that case the solution applies of 
course only to that limited range of m values 
but Ap, (m) 
obtained separately. It 


can be 

that 

A,(@) is equal to the reciprocal value of the 
1 1 


1) for m co 


Ap, (m 
can be shown 


smallest root of the algebraic equation (vide 


p. 941 [6]) : 


L 2 Bo x, t) 0 (11) 


b i 
The reciprocal value of this root is very close to 
the 
u(m) if the corresponding B, is small. 


relative volatility occurring in 
The 
component possessing this relative volatility is 
the 
engineering purposes it is suflicient to use the « 
If a 


more accurate value is considered necessary one 


maximum 


called light key component. For most 


of the light key component for A, ©. 
must solve equation (11). The best way to do 
this is to insert a value of ¢ slightly smaller than 
the reciprocal value of x for the light key com- 
ponent and calculate the corresponding value 


of L,. 
root is then constructed from which the latter 


A graph of L, as a function of the smallest 


can be read for any arbitrarily chosen value 
of L,. 
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Correction for unspecified components 
The absorption factors must now be corrected 
for the presence of unspecified components. The 
been 


complete absorption-factor product — has 


given in equation (5) and the factor C, occurring 
given in 
This factor varies with m but it 
m(AM->.) 


for the unspecified component uw ts 


equation (6). 


can be treated as a constant since ( » a Is 


negligibly small for small values of m. If terms 


containing a negative power of z, are neglected 


and Vy is set equal to D. as a first 


’ 


nt, 
approximation one obtains [4] : 
Ap(M ncaa 1) 


tp, (M+ 1) D, 


(12) 


Then the 
constant of the reference component at the feed 


reciprocal value of the equilibrium 


becomes : 
1, Ap M1) 
/, Apt M) /, 


tray 
I, fp » 1) 
Ky ' V 

. pl. ) 
This expression rmiust hoe equal to the corres- 


ponding expression calculated from the top 


section and in this wa the position of the feed 
located for a given set. of 
,and L,. Instead of A, : 


is used for stripping and top section 


plate ts 
fe 
€é and »"' 


respectively since in case of a variable reflux the 


the notation 


expression analogous to 13 is no longer equal 
to A,' [6}. 

EXAMPLI 
We shall now calculate the absorption factors in 
the stripping section for the following separation : 


Table 1. 


Separation scheme for a sia component system. 





Dist. Kes, x Lib 


Comp. it F 


Ooo SOO) Ton 
oo 
13s 
1M 
Ooo 
OOSTS 


ee woz 
000 
oor 


or 
oie 
M120 
oie 


so) 
1H 
bd 
OS 


(628 i eed 





Product rates in mole sec. Comp. 4 is reference component 


boiling point feed is assumed. 


values of 
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First the curve u(m) has to be drawn which is 
done in Fig. 1 by plotting the lines B, 2" on a 
logarithmic scale as a function of m. 


} 


ulm) i 


"0.155 (age) 


~ 


OMS (235) -—_.---_| i 


‘0.002 (h5e)™ 


. 


r _—>— —_--—-- 


+ -0.16(a.0572) 


—*0.12 (0.302\""———— 
valbetiaaalbadenased 


| | 
10 4 l : | ‘ J 
0 10 


—_—- /7) 





Fig. 1. Characteristic function for stripping section. 

The curve u(m) is then approximated by 
” + : > 

By, %, finite 


of m values i.e. for m 1.2... 10 in this case. 


u* (m) By, ty)" tor a range 

The data pertaining to the substitute com- 
ponents are given in Table 2. 

The MeCabe-Thiele diagram for this system ts 
shown in Fig. 2. The coordinates are the con- 
The operating 

0-0846 which 
is the concentration of component I in’ the 
residue. The ideal trays have been stepped off 
for an operating line corresponding to L,* — 3-313 
and the absorption factor of component 1b on 


centrations for component 1b. 


line intersects the line y = a ata 


a tray is Ltr V,*y since y « is equal to the 
equilibrium constant of the component considered. 
The absorption factor of the reference component 
is then ,, L,*2 V*y and this is also its absorption 
factor in the multicomponent-distillation for 
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Table 2 


Substitute components for stripping and top section 
stripping section. 





Comp. Res. Lib 
1b O-O1L45 1:35 
2h 0-155 O84 

u* (0) 0-1695 





By, u*(0) = 0-0846 ; 
0-202, 


%1p/%p = 1-61 ; u(O) + u*(0) = 


Top Section. 



























































Comp. Dist. “1 
it O215 2-33 
21 0-003 1-035 
u* (0) 0-218 
Hyp /%qye = 225; By, u*(O) = 0-985; u*(0)—u(O) = 0-410. 
900 
a: | 
3 Wy 
& £ y Uj 
a: 
. 
et & 
” 
60F ZA 
F 
r y, 
e UF fo 
3 oof te = 10.0 
20k 
’ f 
0 Pppptittsi disp ss sips tisisssiiil 
0 20 £0 60 80 100 


Mole % Liquid of Comp.b —= 


Fig. 2. Separation of substitute components 1b and 2b in 
stripping section, relative volatility = 1-61, 
which L, = L,* + 0-202 3515 for m = 2, 


%....10. At the first tray however we have an 
absorption factor of the multicomponent system 
equal to (R, R,*) Ap,*(1) = 0-454 Ap,* (1). 
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The absorption factors and absorption factor 
products are given in Table 3. 


Table 3 
Absorption factors and absorption factor products for 





Rp = 400, Ly* = 3-313. 
m Am Ap (m) 
1 0-420 0-420 
: 0-960 0-403 
3 0-977 0-394 
1 1-02 0-403 
i) 1-07 0-430 
6 1-12 0-481 
7 1-16 0-559 
8 1-21 0-677 
9 1-25 0-847 
10 1-28 0-109 





The limiting value of Ap, (m), Ap, (m — 1) for 
m — oc is then obtained as the reciprocal value 
of the smallest root of equation (11). The function 
€ is now calculated from equation (13), thus it is 
successively assumed that the trays M+ 1 = 5, 





















































ae are the feed tray and the curve € as 
a function of M +1 is drawn in Fig. 3. The 
curve is not extended to an M + 1 smaller 
25 oy a 
Sory 
20 “7 
——— 
5 = 
2 
107 — 
*l or Nel 
Fig. 3. Values of £ and » corrected for unspecified com- 


ponents for six components system as a function of M+1 
and N+1, 
for Rp = 40 
for Rp = 10-0. 


than 5 since the residue rate can only become 
negligibly small for an unspecified component if 
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the number of trays is so large that the term 
Ap(M +1) 2," C, 
1 is small. 

Anticipating the discussion of the top section 
it might be mentioned that » in Fig. 3 is the 
reciprocal value of the equilibrium constant of 
the reference component at the feed tray calcu- 
lated from the top section. At the feed tray we 


-M oceurri in e 
a, occurring In equation 
(5) for m 


must have » = € and thus all possible locations 
of the feed tray can be read from Fig. 3. The é 
and » curves corresponding to another reflux 


ratio are also given in Fig. 3. 




















—/.. 


Fig. 4. Liquid- and vapour-rate of component 5 in the 
stripping section. 
Representing VV; ,, = Li mss > B, 
Representing V; ,, = % Sp Lim: 


Liquip Rare or A COMPONENT 
The liquid rate of the lighter substitute component 
can be directly read from Fig. 2. We have 
L, for the liquid leaving the 
The stepwise construction of L,, ,, 


Lin .m "bom 
mth tray. 


consists essentially in alternate application of 
the equation of continuity. 


V.. =L,.., — Bp (14) 


1m 1m 


and the equilibrium condition : 


Vv, = 4,5, L, .. (15) 


mm 


Equation (14) is represented by a straight line 
with a slope one in a diagram in which the 
coordinates are V, ,, and L, ,, and equation (15) 
is represented by a straight line through the 
origin with a slope «,; S,,. 

Since §,, is different for each tray one obtains 
a number of such equilibrium lines each repre- 
senting the relation between V; and L,; on a 
certain tray. Thus it is now possible to obtain 
L; ,., by stepping off from the operating line 
to each of the equilibrium lines in succession. 
In Fig. 4 this procedure has been applied to 
component 5 in Table 1 for a feed tray location 
of M + 1 = 6. The correction for the unspecified 
components according to equations (5) and (12) 
has been applied to the absorption factor products 
in Table 3 and then S,, = Ap(m — 1) Ap(m) 
is calculated. 


The top section 


The equation in the top section equivalent to 
equation (4) for the stripping section is : 


ken 

VSp(n) = 2 2 D,a*" Sp (k) + 
k=0 

"+ Sp(N + 1)L,L,7 2,” 


(N+1)-n (16) 


vv, 
ed | Lio x; 


Lam +1 %u 


The product of stripping factors of the reference 
component is denoted by Sp (n) 


Sp(n) 


and the trays are numbered downward starting 
with n = 0 for the condenser. In the top section 
components 1 to 4 inclusive are specified and 5 


S,S,....8, Sp0)=1 (17) 


and 6 are unspecified. 
The solution of equation (16) is : 


Sp(n) 


in which Sp,(n) is the solution of equation (16) if 
all unspecified components are neglected and C,, 


Sp,(n) +Sp(N +1) | ag Cx mM>1) (18) 


Is : 
L, LY, | (18 
lL, V, 2; a,"* u (n k) 


C,, (for top section) ») 
in which u(n — k) is written for 2, Da" 
One obtains in a similar way as in the stripping 
section the stripping factor of the reference 


150 





VOL, 
3 


1954 








W.R. van Wisk and H. A. C. THtsssen : 





component corrected for the presence of unspeci- 
fied components : 


Sp (N T 1) Cc. a HN+D) 
Sp, (N + 1) B, 4 
2.’ BL, “ 
Ll, 


-n>~(N~D) (20) 


In this equation positive powers of x, have been 
neglected in the expression for Sp(N + 1) C,. 
Equation (20) is entirely analogous to eq. (12). 
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10 4| (6 ) 1: 
0 
_—s 
Fig. 5. Characteristic function for top section. 


The procedure of approximation of u (n) of the 
multicomponent system in the top section by two 
substitute components runs along entirely ana- 
logous lines as that for the stripping section. This 
is done in Fig. 5 and the data have been collected 
in Table 2. The McCabe-Thiele diagram for this 
system is shown in Fig. 6. The stripping factors 
of this binary distillation are the same as those 
of the multicomponent mixture without unspeci- 
fied components from x = 2 upwards to n = 10 
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but a vapour rate V* = V, — u(0) + u* (0) in 
the binary distillation corresponds to V, in the 
multicomponent distillation. We have further 
Spl) = {Rp R,,*}Sp*(1) and for infinite » the 
stripping factor becomes equal to the reciprocal 
of the value of the smallest root of the algebraic 
equation : 

V 2) DO 


at) =0 (21) 


t 





100 


_oe 








60: 


vapour of Comp. t, 


ec 
70 





Mole 








[_— 7» Om 60» §=6©. 80 100 


Mole % Liquid of Coma t => 


Fig. 6. Separation of substitute components 1f and 2¢ in 
top section, relative volatility = 2-25. 


The reciprocal value of the equilibrium constant 
of the reference component at the feed tray K 7 
becomes : 
cas V, SAN) ia: _Sp,(N) ; 

L,Sp(N+1)  L, Spy(N+1) 


(1 — 2,” 3.) (22) 
Ly 

These values have been plotted in Fig. 3, The 
stripping factors and stripping factor products are 
given in Table 4. 
obtaining V;, ,, 
which is similar to the method of obtaining L, ,, 
in the stripping section. These values could 
eventually be used for improving the correction 
for the unspecified components by application of 


A graphical procedure for 
in the top section can be used 
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Table 4 
Stripping factors and stripping factor products for Rp, ~ 4-00 
L,* = 2-53. 

m Sn Sp (n) 
1 0-164 0-164 
2 0-480 O-O787 
3 0-497 O-O391 
4 0-530 O-0207 
5 0-582 O-Ol21 
6 0-654 0-00789 
7 748 0-00590 
S O-833 0-00492 
rT) 0-900 0-00442 

10 O-o1s 0-00404 





the original formulae, equations (6) and (19), 


respectively. 


NOTATION 


The standard system of nomenclature for chemical 


engineering operations (Trans, Amer. Inst. Chem. Engrs. 
1944 40 251) has been followed as far as possible. * Indi- 
cates that the quantity refers to the binary mixture 
approximating a multicomponent mixture. 


- relative volatility of component towards a 


xj 
reference component 
A,,, = absorption factor for component ¢ on tray m 
A,, = absorption factor for reference component on 
tray m 
Ipim) = Ay, Ag... A,,. absorption factor product of 
reference component on tray m 
B = residue rate, moles sec 
BR, = residue rate of component i, moles sec 
6 = subscript referring to bottom section 
C = constant used in correction for unspecified 


components 


D 

D, 

’ 

Kin, 
Lin» Ley 
L,. Ly 
Lim 

m 
M1 
” 

N 
N+1 
Ry, Ry 
Sin 

S, 

Sp (nm) 
“u(m) 
u* (m) 
u(n) 
u®™ (n) 

y os Vy, 
Vas Vy 
V in 
im: “iD 
Yim: 4iD 
7) 
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Summary—The influence of liquid concentration and reflux ratio on plate efficiency has been 
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studied for the system n-heptane-methylcyclohexane using a sieve plate column with 8 trays. 
A definite decrease in plate efficiency is found if the concentration of one component is of the 
order of 0-1 mole per cent. but the decrease in plate efficiency is less than reported for the mixture 


ethanol-water. The system n-heptane-methylcyclohexane appears to follow Raoult’s law better 





Résumé 


than would follow from the existing equilibrium data. 


En utilisant une colonne adiabatique a huit plateaux avec le systéme heptane-normal- 


methyleyclohexane, les auteurs ont étudié Vinfluence de la composition sur lefficacité des 


plateaux. 


Lefficacité diminue lorsque le titre molaire d'un des deux composants s‘abaisse au-dessous 


de 0-1, 
eau-ethanol ou eau-propanol. 


Cette diminution reste pourtant inférieure a celle qui a été signalée pour les systémes 


Le systéme heptane normal-methyl cyclohexane se conforme mieux a la loi de Raoult 


que l'on ne serait amené a le penser d’aprés les données d'équilibre publiées. 


INTRODUCTION 
The question has often been raised as to whether 
plate efficiency in distillation columns depends 
upon the concentration of the components in 
the mixture. 

Keyes and Byman [1] found a sharp increase 
from a value below 5 per cent. to approximately 
85 per cent. in plate efliciency with ethyl-aleohol- 
water mixtures in a bubblecap column, if the 
alcohol concentration increased from 0 to 10 
mole per cent. They found a decrease in efliciency 
near the azeetropic composition. A rapid decrease 
in Murphree plate efliciency near the azeotrope 
of isopropanol-water 
Laxncpon and Keyes [2] (also in a bubblecap 
CHAKVEKILIAN [3] found a 
variation of the H.E.T.P. in a 
column for the carbontetrachlorid- 
benzene. The number of ideal plates varied from 
4 at 10 mole per cent. CCI, to 2 at 80 mole per cent. 
In these investigations, however, the slope of 


has been reported by 


column.) linear 
wetted wall 
system 


the equilibrium curve the liquid viscosity, the 
surface tension, the volumetric liquid velocity 
and the time of contact varied markedly at the 
different compositions. The reported influence 
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of concentration might therefore partly be due 


to other factors. Further a small error in the 
location of the equilibrium curve produces 
inaccuracies in the calculations of the plate 


efliciencies, which become extremely large in the 
neighbourhood of the azeotropic composition. 
Separations of isotopes [4. 5, 6] show that in a 
wetted wall column the H.E.T.P. remains finite 
at extremely low concentrations but no definite 


High 


of plates have been reported for the concentration 


values have been reported. numbers 
of flavours from apple juice although the com- 
ponents involved are very volatile [7]. 

In the 


juice extremely low concentrations of the flavour 


distillation of flavours from fruit 
bearing components, for example, of the order 
of 1 p.p.m. or less occur and a decrease of plate 
efliciency for low concentrations would be very 
important for practical applications. It is for 
that the 


investigated the influence of concentration on 


this reason present authors have 
plate efliciency in a specially constructed column 
in which it has been possible to determine the 
efliciency of the individual plates and to study 


the other factors which might influence plate 
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efliciency separately. The system chosen for the 
experiments was  heptane-methyleyclohexane, 
various reflux ratios and concentrations have been 
used. The properties of this system which is very 
suitable for the present purpose are collected in 


Table L. 


Table 1. Physical properties of the system n-heptane- 


methyleyc oherane. 





Viscosity Slope of the 


equil, curve 


, n-heptane Spec. gravity 


Ho F oO F (abs) ep 


774 1076 
O7S1 


0-688 0-931 





KQuirMENT 


A column was used composed of nine sections, 
four of which had a glass wall; the other five 
were of metal and were provided with a device for 
measuring liquid tlow rate and for sampling. Each 
section except the highest one contained a sieve 
plate. Liquid flow-rate, temperature, composition 
heat the 
measured in a measuring section. The condensor 


and transfer through wall could be 
and reflux divider were mounted directly over the 
top section. The measuring sections were insulated 
by 1-5em of glass blanket, the glass sections 
by a double wall. Every section was surrounded 
by a jacket constructed partly of glass. the tem- 
perature of which could be regulated by means of 
For the 


a set of heating coils. measurements 


of the liquid rates on the trays a new type of 


siphon-flowmeter was developed which combines 
small dimensions (7 em) and a measuring range of 
O-O1-1l ce see 3%. 
minimum 


with an accuracy of about 


designed for 


Sample taps specially 


sample holdup were provided on all measuring 
sections and on two of the glass ones. Additional 
details regarding the tray and column design 
are given in Table Il. A complete description 
of the column will be presented in a subsequent 
paper. 

For 
runs under partial reflux a part flowed by gravity 


Overhead vapour is wholly condensed. 


Concentration and plate efficiency in distillation columns 


Table. 11. 


Summary of plate dimensions. 





Inside diameter of column, em 
Plate spacing, em 
Holes plate 

Hole diameter, em 
Drilling pattern Equilateral triangular 
Distance between hole centres, em 

Weir height, em 


Inside diameter of downpipe, em 





via a pilot valve and a siphon-flowmeter to the 
reboiler. The flowrate of cooling water in the 
the aid of a 
rotameter, and the current in the electric heating 
plate of the reboiler stabilized to within 0-5%. 


condensor was controlled with 


The system distilled was prepared from pure 
grade n-heptane of the Phillips Petroleum Company 
and pure grade methylevclohexane of the British 
Drug Houses. The index of refraction remained 
unchanged after percolation over dry silica gel. 

A continuous record of the vapour tempera- 
ture of all sections was obtained by means of a 
Brown 12-point recorder. 

The 
determined by 
The 
by Prof. vax Hees and coworkers in Delft [8]. 


the 
the 
has 


compositions — of samples were 


means of refractive index. 


refractometer used, been developed 
With this apparatus concentration differences 
of + 0-02 mole % n-heptane could be indicated, 

The refractometer was directly calibrated 
by test samples composed of the pure components 
The refractive indices themselves 
Samples of 0-3 to 0-4 ce 
were drawn; with this quantity no interference 


by weighing. 


were not calculated. 
with the operation of the column was observed. 
The liquid remaining in the sample line was 
blown back into the column just before sampling. 
After steady-state operation was indicated by the 
constancy of temperature and of the composition 
of the samples at the top-plate a complete set of 
measurements was made every two hours. The 
jacket were adjusted 


temperatures every 10 


minutes. 
PROCEDURE 


both at total reflux 
partial reflux with reflux ratios ranging from 


Runs were made and at 
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L, V O-5 to L 
total retlux covered almost the entire composition 


= 0-8. The sets of runs at 


range of the binary system n-heptane-methyl- 
evelohexane. The diffusional driving force 
4* ns, —¥, ranged from 0-085 mole %, at 1 or 
at 99 mole % n-heptane to 1-8 mole % at 50 
mole % n-heptane. 

The superficial vapour velocity for the total 
reflux runs at 50 mole % n-heptane was varied 
in regular increments from 10 to 27 cm. sec. 
The lower limit was determined by the lowest 
distillation rate that permitted stable tray 
operation ; this lower limit or weep-point, was 
carefully established by observing through the 
four glass windows the lowest rate at which 
the trays did not yet start leaking. The upper 
limit was determined by the dimensions of the 
siphon-flowmeters. Entrainment appeared to be 
still negligible. To avoid the effect of a varying 
vapour rate on plate efliciency, the runs at partial 
reflux were made at a single superficial vapour 
velocity of 17 em sec, which is in the centre 
of the range of stable operation. The distillation 
rate of the other total reflux runs was also main- 
tained on 16-5-19-3 cm sec. 

The vapour velocities in the holes amounted ; 
to 238 cm sec for the minimum rate to 642 cm_ sec 
for the maximum rate and to 405 cm. sec for the 
runs at partial reflux. That for velocities in the 
holes ranging from 238 to 500 cm sec bubble 
formation rather than tube flow occurs, can be 
shown from hydraulic considerations. 

In order to eliminate the influence of possible 
traces of volatile impurities or water a small 
amount of distillate was withdrawn slowly as 
soon as steady state operation was attained. 
After this a sample was taken every 15 minutes. 
In each run except, 12, 13 and 19 to 27 three 
or four liquid samples of the trays 1-3-5-6-7-8 
and the overhead vapour were taken, and three 
vapour samples of tray 6. For the runs 12, 13 
and 19 to 27 it was unnecessary to analyse 
each tray on account of the very small difference 
in molal composition between successive trays. 
In these runs only 7 samples of each the lowest 
and top section were taken. At the end of each 
run the liquid flowrates of tray 1-3-5-7 and the 
top section were determined. 
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CALCULATIONS 
The Murphree plate efliciency which depends 
upon the intensity and type of mixing of the 
liquor on and of the vapour below the plate 
as well as on the term mV L was related by W. K. 
Lewis [9] to the Murphree point eflicieney which 
is essentially independent of such factors. It 
must, however, be remembered that neither 
the plate efliciency nor the local or point efficiency 
is a valid representation of plate performance. 
The real values of point efliciencies lie in a range 
between those of the plate and point efficiencies. 
Since for a bubble plate 17}in. LD. the real 
local efliciency appeared to lie just between 
the plate —- and point efficiency [10] and taking 
into account that the caps impete the mixing, 
it may be expected that the plate- and point 
efliciencies in this study are almost indentiecal. 
The overall point efliciency on a gas basis is 


defined by the following equation : 
' | ™ 
Ek Yn, Yn 1 (1) 


For the system  n-heptane-methyleyclohexane 
the compositions of the vapour entering and 
leaving a plate are nearly equal especially at 
the low- and high composition range. Under such 
conditions the subtraction of two numbers of 
nearly the same value is likely to give questionable 
results. In order to avoid this, Murphree point 
elliciencies are calculated from the number of 
overall gas phase transfer units of the entire 
column. As for the runs at total reflux the integra- 
tion is carried out over a very narrow concen- 
tration range (slope of the equilibrium curve, 
diffusional vapour driving force and physical 
properties are hardly varying) this procedure 
may be substituted for a determination of 
individual plate efliciencies. The number of overall 
gas transfer units is defined as follows : 


Yn Yn-1 
\ dy sf (2) 
or jJy—v idss-¥ 
ms wn 


In the interval y, , to y, the integral 
= y —y* 


varies from to Its average 
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Table 111. 


Concentration and plate efficiency in distillation columns 


Summary of the total reflux runs. 





Run number 

Molal liquor-vapour ratio 
Superficial vapour velocity em /sec. 16-8 
Overpressure in the reboiler. cm of water 10-0 


Liquor composition, (mole °, n-heptane) 
of tray number : 


— 
~ 


yo go 
P& as 
wt S&S we he te 
wo 


Molal vapour composition top tray, °, 


~- 
_— 
~ 


0 


Diffusional vapour driving force mole 
(Yn+1 Y,) 


Per cent Murphree efliciency 


calculated according to eq. 5 


plate 


Accuracy of Murphree plate efficiency 
(per cent) 


13-8 


95-01 48-36 
51-36 
54°36 
55°85 
57-32 
58-82 


98-52 
98-60 
98-7: 

YS-S1 
98-90 
98-96 
96-64 


99-02 60-27 


0-085 0-267 0-276 1:77 





+ 


if the 
4, 1 y*.4 1 4, 


value is 


. * 
2 Ls n 
variation is considered as linear in 


this small 


interval. Inserting this average value in the 


integral in eq. 2 gives: 


N OCT ~~ 


0-5 oy, (3 
a : Y, Y,, De Y,, : * y,, . ) 


Y n-l 


substitution of equation (1) in (3) gives 


, - a ’ ’ Yy, 4», 1 
Noor 05.2 ((Eoa, + Fou... x ) (4) 


. Yet Yn 
If it is assumed that E,, is independent of the 
term (y*, — y,,-;) and the slope of the equilibrium 
curve and slope of the operating line are both 
constant over the range of compositions con- 
sidered, the relation between E,,, and Nog, may 
now be written as follows : 


Nn. 
0G, 5 
0G 0-5 .n(1 + C) (9) 


in which C = 9» — 91 
Yn+i - Yn 


can be determined 


graphically. A plot of equation (5) is given in 
Fig. 1, where E,,, is plotted against Nog, with 
the term C as parameter. Obviously the absolute 


value of Ng, is not representative for the 














Fig. 1. Relationship between Ep,, overall gas plate 
efficiency, and Noe, number of overall gas-phase transfer 
Y, — Y, 
units of the entire column, with C = —"—— ot 
n+l ~~ Yy, 
parameter according to eq. (5). 
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Table 111, 





Concentration and plate efficiency in distillation columns 


Summary of the total reflua runs (cont.) 










Run number 19 20 





Molal liquor-vapour ratio 10 10 


Superficial vapour velocity em /sec. 16-8 16-5 





Overpressure in the reboiler. em of water 10-0 10-0 





Liquor composition, (mole °%, n-heptane) 
of tray number : 


1 
”. 
” 
- 
~ 


8 


Molal vapour composition top tray, %, 2-96 5-86 


Diffusional vapour driving force mole %, 


(Ynat 4) O18 O34 


Per cent Murphree plate efficiency 
calculated according to eq. 5 70-5 76-5 


Accuracy of Murphree plate efficiency 
(per cent) 2-2 1-2 





1-99 3-82 






21 


oe 
-- 








1 I) I) 1 1) 10 0 







16-5 16-5 16-5 1-33 19-38 18-7 le-7 





1 li) 1-0 





4-45 +11 7-28 V7 -77 V6-89 95-41 93-18 





6-92 74 1108 9855 ORAM 9T-OG | «95-65 
O-41 O45 0-62 O123 O169 O246 O36 
770 s0O-0 785 75-4 TO4 S17 83-2 
1) 1) eu 3 2-7 17 10 





Murphree efficiency since C decreases for increasing 
values of VL, therefore Nour 
constant values of E,,. Equation 5 is rigorously 
valid for the runs at total reflux in this study, 


since (y*, — y,-;) as well as m are constant. 


decreases at 


The vapour concentrations on individual 
trays can be calculated from eq. (5). In those 
sases in which these concentrations were measured, 
excellent agreement was found with the experi- 
mental values. The Murphree point efliciencies 
of the runs at partial reflux were calculated both 
from equation (5) and the individual tray com- 
positions. Vapour compositions in the partial 
reflux runs could be directly calculated from 
liquid samples and flowrates. The discrepancy 
of the molal liquid flowrates between the trays 
1, 3, 5, 7, and the topsection amounted at most 


to 3% which induces a maximum error of 15% 


in the position of the operating line. In all 
runs the column could be kept adiabatic to within 
05%. 

The vapour liquid equilibrium data of 
Bromitey and QvuiGGLE [11] were tested by a 
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modified Gibbs-Duhem equation of Rerpiicn 
[12] and appeared to be highly inconsistent 
in the concentration range from 80-100 mole 
per cent. n-heptane. On our request the equili- 
brium data were redetermined by the “* Central 
Institute of Physical and Chemical Constants, 
Netherlands.” The new data fitted Raoult’s law 
to within the accuracy of the experimental error 
(0-1 mole ° 
with the experiments of Crorzex, Haase and 


» n-heptane), which is in accordance 
SiG [13] who could not detect any deviation from 
ideality in the liquid state above 40°C. In view 
of this it was considered that Raoult’s law 
could be applied with confidence. The vapour 
pressures of n-heptane and methyleyclohexane 
were obtained from Circular No. C 461 of the 
National Bueau of Standards. 


Resuurs and Discussion 


Total reflux, Table HIT summarizes the value 
of plate efficiency along with the compositions 
of the trays, superficial vapour velocity and 
overpressure in the reboiler, Liquid compositions 
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In order to examine whether this drop in 
etliciency near the extremes of the concentration 
range is not caused by a small deviation of the 
system from ideality, the following testing method 


was introduced. 














%o Murphree plate efficiency —e 


Relative volati/aty —e 


Superficial vapour velocity cm/sec —e 
Fig. 2. Effect of vapour velocity on plate efficiency 
at total reflux Y*,,., — Y, = 180 mole per cent. 


in the table are an average of three analyses. The 
effect of column throughout upon the local i ee pe 
Murphree efticiency is shown to be small, ef. 1954 
- , “ates ey Mole % n heptane in liquid —= 
Fig. 2 which is in general agreement with the ; - 
. : . Fig. 4. Experimental relative volatilities if effciencies 
literature on this subject. Bix = 
M ; aE x ar are independent of the liquid composition. 
Murphree efliciencies at a constant superificial A x ideal 
vapour velocity of 16-5-19-3ft sec are plotted @ x experimental 
in Fig. 3 as a function of molal liquid composition. 
It is apparent that the efliciency varies markedly Assuming that Raoult’s law holds in’ the 
with composition in the ranges of low and of neighbourhood of 50 mole % the values of Murphree 
high concentrations. This variation can not be — efficiencies in this range were calculated. From 
attributed to a lack of precision in the calculation these a new equilibrium curve was calculated, 
of the Murphree efficiency. (See the list of assuming that the efliciency remains constant 
accuracies in Table III). in the entire concentration range. The relative 


"I 





7; 
log z 





& 





3 





%o Murphree plate efficiency —e 





605 





Mole Yo n heptane in liquid —e 
Mole Yo n heptane in liquid —= 


Fig. 3. Effects of liquid concentration on plate efficiency 
at total reflux. Superficial vapour velocity 165- Fig. 5. Thermodynamic testing graph of the new 
19-3 cm/sec. equilibrium curve. 
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volatilities of this new curve are plotted along 
with the ideal volatilities against the liquid 
composition in Fig. 4. Repiicn’s modified 
Gibbs-Duhem equation was now applied to 
test the consistency of this calculated equilibrium 
curve. In Fig. 5 the logarithm of the ratio of the 
activity coeflicients is plotted against the liquid 


composition. The values of log “1 are obtained 
Y2 
from the following equation : 
> 
"1 P, 
log ~ log x log P 
Y2 2 


As the curve does not conform to | log "dr =0 

. Ya 

0 
it must be concluded that the variation of the 
Murphree efficiency with concentration can not 
be attributed to either inaccuracies or wrong 
equilibrium curve and that an influence of the 
diffusional driving force upon efficiency exists. 
A similar effect for the system ethanol-water 
(Keyes and Bymay, [1]) is represented in 
Fig. 6 where E,,, is plotted against liquid composi- 
tion with vapour rate as parameter. The influence 
of the diffusional driving force appears to decrease 
at increasing vapour rate. Owing to the limited 
accuracy of sample analysis and varying physical 
properties with concentration no further con- 
clusions will be drawn. 


Yo Murphree plate efficiency 





Mole °/o ethanol in liquid —.» 


Fig. 6. Effect of superficial vapour velocity and liquid 
composition on plate efficiency for the system ethanol- 
water [1]. 
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Harcuer and Sace [14] in an investigation 
of the influence of natural and forced convection 
upon the transfer of n-butane to and from hydro- 
carbon liquids found for natural convection 
conditions a fairly regular decrease of the mass- 
transfer coefficient for the liquid film as the 
concentration of the n-butane in the liquid 
phase approached that at the interphase, while 
this effect appeared to be negligible with forced 
convection, which is in accordance with Fig. 6. 

Partial reflux, In the runs with partial 
reflux the distillate was returned to the reboiler 
in order to obtain a stationary state. In Fig. 7 
the analysis of the individual liquid samples of 
partial reflux run 17 are plotted against their 


corresponding tray numbers. In the same graph 


TE a 
| 
| 


re) 
% 
—= 
| 
— 











Mole %o n heptane in liquid —e 








Tray number —e 
Fig. 7. ‘Tray compositions for a typical partial reflux run. 
@ By liquid analysis 
— Calculated curve I for Egg = constant. 
Calculated curve II for Egg = f(¥*n.y — Yn) 


curve (1) represents the liquid compositions caleu- 
lated according to the BAKER-STOCKHARDT pro- 
cedure [16] so constant plate efficiencies over all 
the trays and a curve (II) representing liquid 
compositions giving plate efliciencies which are a 
same function of the diffusional vapour driving 
force as found from the total reflux runs at the 
same vapour rate. From a separate study it 
appeared that the samples of tray 3 and 5 might 
have contained some condensed vapour (too high 
concentration of n-heptane). It is therefore 
concluded that curve IL must be preferred. A 
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downward trend of efliciency with decreasing 
reflux ratio is indicated in Fig. 8. 

It has been stated that the plate efliciency 
still increases at increasing values of LV above 


unity in the bottom section (15). In that case 


, 





3 
= 
2 80 
= 
aq 
v 
a 
= 70 
a 
> 
50 + 
/ eee eee ow re bi 4 sp ss 9J 
C 06 07 08 09 10 


Reflux rato L/V —e 
Fig. 8. Influence of reflux ratio on plate eflicieney. Super- 


ficial vapour velocity 17 em sec. Average concentration 
of n-heptane in liquid 350 mole per cent 


the relationship between Ay* and FE) cannot be 
explained only by the effect of a simultaneous 


decrease or increase of the diffusional driving force. 
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It is therefore concluded that the plate efficiency 
is dependent on the term mL V, according to the 


relation between the overall-and individual 
3 : 1 1 mVOL 
number of transfer units _. _ + . as 
Nog Ng Ni 


well as on the diffusional driving force. 
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Organisation 


NOTATION 
Ep, — overall point efficiency, based on gas concentra- 
tions. 
Ng = number of gas-film transfer units. 
N;, = number of liquid-film transfer units. 
NOG, overall number of transfer units, based on gas 


concentrations. 
L_ = molar liquid flow (g mol/sec). 
J molar vapour flow = (g mol /sec). 
m = slope of equilibrium curve dy* dz. 
" number of plates. 
P Vapour | ressure of pure comy onent 


x = relative volatility. 

y = activity coeff. 

y = mole per cent. in vapour. 

t mole per cent. in liquid. 
y* equilibrium value of y corresponding to x. 
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Extraction of uranyl nitrate in a disc column 
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Summary —The dise column designed by SrerHeNs and Morris [14] for gas absorption studies 
has been modified for countercurrent liquid-liquid extraction. The effect of liquor rates on the 
overall mass transfer coefficients and heights of overall transfer units has been studied for the 
continuous, countercurrent extraction of an aqueous solution of pure uranyl nitrate. The feed 
solution contained 0-348 Ib. UOg(NO 3) per Ib. of solution and the extractants were uranium-free 
diethylene glycol di-n-butyl ether (dibutyl carbitol) and methyl iso-butyl ketone, initially con- 





taining 0-049 Ib. UOg(NOg)y per Ib. of solution. 


each case. 


The solvent was the “ dispersed’ phase in 


The data may be correlated, for the range of liquor rates studied, by the equations : 


(a) Dibutyl Carbitol Extraction : 
Kp 1-20 R,,” 24 E,,° 76 
Kp = 085 R924 B07 
0-76 


R 
(HTU oR vs ” 
E,, 





R “O24 
(MTU jog 0-7 = 
E,, 


(b) Methyl Iso-butyl Ketone Extraction : 


Kp = 421 R,,°* £0 


Kp = 2-91 R,,°* E,,°* 


R 0-64 
(HTU) oR = 1-95 = 
Ey, 

~) 86 
(ATU og = ae | _ 
“mm 


Résumé — La colonne a disque congue par SterueNs et Morris pour l'étude de absorption 
gazeuse a été modifiée pour servir 4 lextraction liquide-liquide a contre-courant. 
L’influence des vitesses de la liqueur sur le coefficient de transfert massique principal et 
sur les hauteurs dunités de transfert principales a été étudiée pour une extraction continue a 
contre-courant d'une solution aqueuse de nitrate duranyle. pur. La solution d’alimentation 
contenait 0,384 Ib, de UO, (NOx), par Ib. de solution et les extracteurs étaient le diéthyléne glycol 
di-n-butyl éther (dibutyl-carbitol) sans uranium et le méthyl iso-butyl cétone contenant 
initialement 0,049 Ib. de UO, (NOs) par lb. de solution, Le solvant était dans chaque cas la phase 


* dispersée.” 


INTRODUCTION 
Continuous, countercurrent — liquid - liquid 
extraction investigations are usually performed 
in wetted-wall, spray or packed towers. With 
the former apparatus the area for mass transfer 
is known and if a glass tower is used only the 
aqueous phase will flow as the *“ dispersed ” 
phase around the column wall. With the latter 
two columns either phase may be dispersed, but 
the area for mass transfer is unknown and 
necessitates the use of * volumetric mass transfer 


coefficients ” or “extraction capacity coefficients,” 
Ka, to correlate the data. 

The dise column has proved of value for 
correlating mass transfer coeflicients with liquor 
rates in gas absorption studies, as shown by 
SrerueNns and Morris [14] and Roper [10, 11]. 
This column offers an advantage over packed or 
spray towers because of the known area for mass 
transfer and is particularly suited for gas-liquid 
system investigations. Experimental data have 
not been published on the use of the dise column for 


* Present address : Chemical Engineering Division, A.E.R.E., Harwell, Berkshire. 


161 















R. K. Warner : Extraction of uranyl nitrate in a dise column 


liquid-liquid system studies. Satisfactory operation 
of the column for such work depends primarily on 
the choice of the disc material. The appropriate 
phase must flow as a uniform film over the complete 
area of the packing element, with essentially no 
“dispersed” phase droplets detaching themselves 
from the dises. With suitable disc materials, 
either the solvent or the aqueous phase should 
tlow satisfactorily over the packing element and 
thus allow flexible operation of the column. 
When this work was begun only two papers 
were available which reported data on the 
performance of extraction columns for the 
extraction of inorganic compounds by organic 
solvents. The extraction of ferric chlomde from 
an aqueous hydrochloric acid solution by di-tse- 
propyl ether in a spray column, when the solvent 
was dispersed, has been studied by GrankopLis 
and Hixon [2]. With this system (H#TU),, 
decreases with increasing ether flow rate and 
tends to a minimum of approximately 0-7ft at 
high solvent rates and ts independent of the 
ferric chloride concentration. The separation of 
cobaltous chloride from nickel chloride by 
extraction of an aqueous hvdrochlone acid 
solution of these salts with octanol-2  (capry! 
alcohol) in a spray tower, with the solvent phase 
dispersed, has been investigated by KyLanprer 
and Garwin [5]. Average cobalt : nickel ratios 
of approximately 200; 1 in the extract were 
obtained from a feed containing about 1-3; 1 
ratio of these salts. The extraction capacity 
coefficients K, a and Ky a increase with increasing 
solvent phase rate and are independent of the 
cobalt ; nickel ratio in the aqueous phase. For 
an essentially constant aqueous rate (TU), 
decreases slightly within the range 6-4ft. with 
increasing solvent rate, whereas (HTU),, de- 
creases rapidly under these conditions. It is 
reported that the overall H7TU values for the 
extraction of this inorganic solute are of the 
same order as those obtained for the extraction 
of organic solutes in similar equipment. 
Merpocn and Prarr [9] recently published 
data on the extraction of aqueous uranyl! nitrate 
solution by dibutyl carbitol and methyl! iso-butyl 
ketone and for urany! nitrate transfer from dibuty| 
carbitol solution into water. A wetted-wall 


column was used and a range of concentrations 
were studied. It was found that the individual 
film resistances are approximately equal, while 
the interfacial resistance. due to the solvation 
reaction, is 5-20°, of the overall for dibutyl! 
carbitol and 10-35°,, for methyl! iso-butyvl ketone. 
The aqueous film mass transfer coetlicients vary 
inversely as the concentration and satisfactory 
correlation with organic solute transfer data was 
obtained for both tilm coeflicients, 

The present work was undertaken to determine 
the characteristics of the disc column for liquid- 
liquid extraction and particularly to provide data 
on the extraction of an inorganic salt from 
aqueous solution. Current interest in the extrac- 
tion of aqueous solutions of inorganic compounds 
centres largely around the extraction of cerium 
(IV) [19, 22], thorium [1, 7, 12, 15, 16] and uranyl! 
nitrates [3. 6, 8, 13) by oxygenated solvents, 
Extraction of uranyl nitrate is commercially 
unportant because of the case of extraction, the 
extremely high selectivity of certain solvents and 
consequent purity of the extracted material. This 
salt may be extracted from its aqueous solution, 
with or without the addition of metallic nitrates 
as “ salting-out ” agents, by a wide range of 
oxygenated solvents [4, 8, 20]. Extraction of an 
aqueous solution of pure uranyl! nitrate containing 
0-348 Ib. VO,(NO,), per Ib. of solution by dibutyl 
carbitol (D.B.C.) was chosen for the initial work. 
D.B.C. was selected because it possesses suitable 
extraction properties and high selectivity ; it is 
readily recovered; the distribution data are 
available and its low volatility facilitates opera- 
Furthermore, D.B.C. is 
comparatively viscous and with a dise material 


tion of the column. 


readily wet by the solvent, it was considered 
uniform solvent flow over the packing element 
would be maintained. To study the effect of a 
different solvent type, further runs using methyl 
iso-butwvl ketone (M.I.B.K.) and essentially the 
same inlet aqueous solution were performed. An 
inlet solvent concentration of 0-049-Ib. VOANO,), 
per lb. of solution was adopted for these runs. This 
was used to avoid the highly curved lower portion 
of the equilibrium line and therefore allow the 
possibility of using the logarithmic mean driving 
force in the calculations. 
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EXPERIMENTAL 
(a) Apparatus 
The arrangement of the apparatus is shown in 
Fig. 1. 
up the dise assembly as the 


The solvent phase was required to flow 
* dispersed ”’ phase 
for this study. Therefore the dise material must 
be preferentially wet by the solvent if a uniform 
film is to be maintained. Carbon discs were chosen 
after preliminary experiments with wood and 


carbon discs showed that with the latter. uniform 


wetting and solvent flow were obtained. Further 


desirable properties of the hard carbon dises are 
chemical stability and low porosity. 
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Fig. 1. General arrangement of apparatus. 


To provide a greater area for mass transfer 
than has been used previously in dise column 
studies, 79 carbon discs, each 0-0484ft. diameter 
and 0-0148ft. used, The discs 


thickness were 
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were threaded edgeways on a thin glass rod and 
maintained at right angles by means of a methyl! 
methacrylate-ethylene dichloride cement. The 
dise assembly was positioned centrally in a lin. 
1.D. 
and the column was mounted vertically. A 
the 


The lines to the column were }in. LD. 


glass tube fitted with suitable connections 


detailed diagram of column is shown in 


Fig. 2. 
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Fig. 2. The dise column for liquid-liquid extraction. 


glass tubing with polythene connections for lines 
containing solvent and neoprene connections for 
lines containing the aqueous phase. 

To provide suflicient head for solvent flow, a 
small air blower, fitted with a by-pass, was 
connected to the solvent feed tank. The aqueous 
uranyl nitrate solution was passed through a glass 


1638 





R. K. Warner : Extraction of uranyl nitrate in a disc column 


capillary tube fitted with a manometer to indicate 
and allow control of the flow-rate. The solvent 
feed was metered through a rotameter. 


(b) Materials 


Pure uranyl nitrate hexahydrate was dissolved 
in distilled water and the concentration adjusted 
to give a feed solution containing 0-348 Ib. 
UO ANO,), per Ib. of solution for the D.B.C, 
extractions and containing  0-346-0-347 Ib. 
UO ANO,), per Ib. of solution for the M.LB.K. 
For the latter 
feed solution was saturated with solvent before 


extractions. runs the aqueous 


admission to the feed tank. Before cach series 
of runs the feed solution was well mixed and the 
UO,NO,), concentration determined by evapora- 
tion of an aliquot and ignition to U,O,. Batches 
of recovered uranyl nitrate solution were equili- 
brated with the appropriate solvent and the 
solute concentration in each phase determined. 
The in the 
solution was adjusted by either U,O, or HNO, 


uranium : nitrate ratio recovered 


addition to give satisfactory agreement with the 


distribution data in Table 1. 


Table 1. Distribution of uranyl nitrate between dibutyl 
carbitol-water and methyl iso-butyl ketone-water at 20°C. 





Concentration*® 
Solvent 
Aqueous phase Solvent phase 
a, 4, 


1. Dibutyl 0-476 
Carbitol OAT 
OSS 

O-281 

0-246 

O-178 


O-396 
355 
253 
122 
O-O74 
Oo 


0-460 
0-418 
0-355 
O314 
0-300 
0-270 
243 
0-206 
0-149 


O-336 
O- 285 
0-182 
O11S 
0-096 
O-O62 
O-O36 
OO1T 
O-005 


2. Methyl 
Iso-butyl 
Ketone 





* Ib. UOQ(NOg)y per Ib. of solution. 


D.B.C. was supplied by Carbide and Carbon 
Chemical Corporation. The material was used as 
received, except for with ferrous 
sulphate to remove peroxide initially present, 
The solvent 


treatment 


followed by washing and drying. 
conformed to the following specification : 


Specific gravity (20 20°C), 0.883 
Colour, water white 

Solvent content, 98°. min. 

Acidity, 0-01", (wow) as CH,COOH 
Suspended matter, nil 


M.L.B.K. was supplied by Shell Co. (Aust.) and 
conformed to the following specification : 


Specific gravity (20 20°C), 0-801 
Solvent content. 99°. min. 

Acidity, 0.01%, (wow) as CH,COOH 
Distillation range, 114-117 C. 


M.1.B.K. was stirred with concentrated aqueous 
uranyl! nitrate solution until the specific gravity 
of the solvent phase was approximately 0-840, 
The M.LB.K. phase was separated, analysed 
gravimetrically and the concentration adjusted 
to give 0.049 Ib, UO,(NO,), per Ib. of solution. 


(c) Procedure 
The flowmeters were checked before cach series 
of runs to ensure satisfactory operation. 

To begin a series of runs, the column was 
tilled 
stand one hour to give complete wetting of the 


evacuated, with solvent and allowed to 


discs and removal of adhering air bubbles. With 


the solvent flowing at the desired rate, the 


aqueous uranyl nitrate solution entered at the 
top of the column. The solvent surrounding the 


discs was displaced and when the column was 
full of aqueous phase the bottom raflinate line 
was opened, Both liquid rates were adjusted to 
the desired values and the liquid interface main- 
tained inside the arm above the top dise by 
adjusting the height of the raflinate overtlow 
weir. Air and solvent trapped at the top of the 
column were removed through the vent. 

Steady conditions were maintained until the 
volume of raflinate phase collected was approxi- 
mately two litres (three times the volume of 


aqueous phase held in the column). The flow 
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rates of both the outlet raflinate and extract 
phases were then obtained by measurement and 
the specific gravity of each phase determined with 
a Westphal Balance. The flowmeter readings 
were recorded, the temperature of the extract 
and raflinate phases noted and samples taken for 
analysis. Duplicate results were obtained for the 
same further 4-500 ml. of 
raflinate and for certain 
triplicate results were obtained after another 
2-300 ml. of raflinate had collected. The flow 
rates were altered and the procedure repeated. 

Since only a limited quantity of uranyl! nitrate 
was available, only three consecutive series of 


flow rates after a 


had _ collected. runs, 


runs could be made, and it was 
recover the salt from the extract and raftinate 
phases. The D.B.C. extract was filtered to remove 
a light yellow solid and then batch extracted with 
distilled water until free of uranyl nitrate. The 
solvent was dried with MgSO, and returned to 
the feed tank. The from 
treatment of the extract was combined with the 


raflinate, twice extracted with benzene to remove 


necessary to 


aqueous solution 


dissolved solvent and evaporated to the required 
concentration. The M.I.B.K. extract was stirred 
with water to reduce its uranium content to the 
desired figure and the raflinate evaporated as 
above. The distribution properties of recovered 
uranyl nitrate solution were checked and when 
satisfactory the solution was returned to the feed 
tank. 

Aliquots of the ratlinate and extract phases 
were evaporated to dryness and the urany! nitrate 
concentration determined by ignition at 950°C to 
U,0,. D.B.C. extract samples could not be left 
overnight precipitation of small 
amounts of a light vellow solid. 


because of 


CALCULATIONS 


For a small section of the column containing an 
area of dA, the transfer 
coeflicients are defined by Equation 1 : 


dN, = Kx (a Ke (y, 


A material balance for this element of the column 


dises overall mass 


v)dA y)dA. (1) 


gives : 


dN, = WI, dx = WI, dy. (2) 


Extraction of uranyl nitrate in a dise column 


These expressions assume that strict counter- 
current flow occurs. Although circulatory currents 
were observed in the aqueous phase, it was not 
possible to correct for this effect, and the above 
expressions were used without modification. 
From equations (1) and (2), considering the 
relationships for the extract phase, and _ inte- 
grating over the length of the column, assuming 
Kg, is constant and using a mean value of Tp. 


y 


E (mean) dy 


Ky (4. —y) 


“ y 


The assumption that the overall mass transfer 
coellicient is constant is not strictly correct, since 
the distribution varies with 
centration in However, 


coetlicient con- 


these systems. this 
assumption is used in practice, regardless of the 


curvature of the equilibrium line [9, 17]. 


Now by definition, 


soma d 
(NTO los (y : y)” 


and 
Z 


(HTU)o, (NTU jog 


dy 
(Ye 


u 


1 


Therefore, from Equations (3) and (5), 


(HTU) op 


Ve (mean) (6) 
Ky 


Similarly, from the relationhips for the raftinate 
phase, 
; Z 
(HTU op a 
(NTU Jor 


I R (mean) 


(HTU)on = 
R 
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The overall mass transfer coeflicients A, and Ay, 
were calculated from the experimental figures by 
the following relationships : 

N 


Ky : ; (Q) 
. Al AY om 


: N 
K, = 4 (10) 
A( AX) oy», 
where N, is based on the uranyl nitrate trans- 


ferred to the solvent phase, and (Ay), and 


(Az),,, are defined by : 


tm 


(Ye y,) 
( Ay) (2 - a. 
Ym Me (NTU ox 


(t r,) 


(NTU on 
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Fig. 3. The equilibrium curve for the urany! nitrate-dibuty] 
carbitol-water system and operating line for Run No. 8b in 
D.B.C, series. 


Fig. 3 shows the distribution curve for the 
uranyl! nitrate-dibuty! carbitol-water system and 
the operating line for Run No. 8b in the D.B.C. 
series. The concentration driving forces given in 
equations (11) and (12) were obtained by graphical 
integration for the D.B.C. series. Fig. 4 shows 
the equilibrium curve for the uranyl! nitrate- 


methyl iso-butyl ketone-water system and the 
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operating line for Run No. 6b in the M.LBLK. 
Despite the use of a finite uranyl nitrate 
the M.L.B.K., the 
logarithmic mean driving force was only applicable 
in the case of Ay and equation (12) had to be 


series. 


concentration — in entering 


used to determine (Ax) 


tm* 





oe " 


9 ? 
~ - 


LB URANYL NITRATE PER LB SOLVENT PHASE 
° 
2 











2 o3 o« 
iB URANYL NITRATE PER LB AQUEOUS PHASE 
Fig. 4. 
methyl iso-butyl ketone-water system and operating line 
for Run No, 6b in M.I.B.K, series. 


The equilibrium curve for the uranyl nitrate- 


The uranyl nitrate transfer, N, lb. hr.. was 
based on the solvent phase throughout and _ is 
more accurately known for the D.B.C, series than 
for the M.L.B.K. series. N, so calculated differs 
from the solute transfer based on the raflinate 
phase because of small fluctuations in the inlet 
aqueous phase rate and also because of the 
analytical errors associated with the small con- 
centration changes in the raflinate phase. Further- 
more, mutual solubility factors alter the raflinate 
concentration and the term (a, — #,) is not a true 
indication of the solute transfer. 

The evaluation of (H#TU),, by equation (7) 
is therefore unreliable and gives scattered results 
markedly different from those calculated from 
However, (HTU),» calculated 


from either equations (5) or (6) are in close 


equation (8). 


agreement and would be identical if the extract 
rate remained constant throughout the column. 
Therefore the (HTU), values in Tables 4 and 5 
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Table 2. Experimental data — Dibutyl carbitol extraction. 








Lulet Inlet 
Run Raffin. Extract aqueous solvent Raffin- Extract 


No. ate conc.* rate rate ate rate 

cone.* Y ky E, R, ky 
zy lb/hr. Ibohr. lb hr. tbo hr. 

la 0-331 0-02 3-18 0-58 3-26 O-70 
Ib 0-332 0-092 3-16 0-58 3-10 0-69 
le 0-331 0-093 3-20 O57 3-18 o-70 
2a 0-324 0-073 3-20 1-59 3-24 oo 
2b 0-324 0-073 3-20 1-60 3-26 2-06 
2e 0-324 0-072 S15 1-61 3-20 205 
3a 0312 0-063 +18 3-00 3-00 3-38 
Bb 0-311 0-064 3-18 3-00 3-18 345 
Be O-311 0-064 $18 3-00 3-11 349 
4a 0-290 0-054 3-18 5-06 3-02 546 
4b O 289 0-055 3-18 5-06 3-05 56 


5a 0-344 0-127 13-90 0: 
5b 0-344 0-127 13-95 0% 
6a 0-341 0-104 13-90 1 
6b O341 0-104 13-95 1: 13-00 1-83 
7a 0-336 80-092 13-98 14-48 3-42 
Tb O335 0-091 13-95 3-00 13-25 3435 
Sa 0-329 0-090 13-90 5-06 13-68 1-98 
Sb O328 0-090 13-02 506 | 1381 501 


13-84 O-t4 
13-04 O-t4 
13-85 1-83 


tw wt 
as 


— 
~ 
= 


8c 0-328 «80-090 13-900 5-06 13-57 5-00 
Va 0-340) «0-110 540 0-57 5-67 0-68 
ob 0-340) =60- 108 5M O59 5-78 O72 
10a) = OB29(O-08 5M 1-59 5-33 1-90 
10b 0-331 0-081 5-30 1-59 555 1-97 
lla 0-317 O-O72 5-20 3-00 07 $44 
lib) «60-319 0-073 5-20 3-00 ov 3-48 
lle 0-319 0-073 5-20 3-00 $95 3-42 
12a) «600-305 0-069 BAS 5-06 4-62 +31 
12b0 «60-308 (O-068 5-20 5-06 4-70 5-32 





* Ib. VO(NOg)y per Ib. of solution. 


were calculated from equations (6) and (8) rather 


than from equations (5) and (7). 


ResULTs 
The experimental data are shown in Tables 1, 
2 and 3. The dise column dimensions are given in 
Table 6 and the calculated data for the D.B.C. 
and M.I.B.K. runs are contained in Tables 4 and 
5 respectively. 


The temperature of the system varied between 
178°C and 22-3°C for runs with D.B.C. and 


between 22-6°C and 27-5°C with M.LB.K. Earlier 
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Table 3. Experimental data Methyl iso-butyl ketone 


extraction .** 





Inlet Inlet 
Run Raffin- Extract aqueous solvent Raffin- | Extract 


No. ate cone.* rate rate ate rate 
conc*. YW» Ry Ey R, Ey 
ry Ib hr. lb/hr. Ub. hr. Ub. hr. 
la 0-327 0-149 3-80 O-s1 3-76 O89 
Ib 0-328 0-150 3-80 O-s1 3-78 0-89 
le O-327 0-150 3-82 O-s1 3-74 O-88 
2a 0-339 0-162 11-74 O81 11-68 0-383 


2b 0-339 0-162 11-83 O-s1 11-87 O-87 


2c 0-339 O-161 12-02 Os! 11-99 0-85 
ba O3138 0127 381 2:19 3:59 2-42 
3b O313 0-126 B85 2:19 3-77 2-45 
tu O-331 0-150 12-10 2-19 12-01 2-54 
> O-331 O-151 12-49 2-19 12-07 2-58 
+e O-331 0-150 12-12 2-19 12-00 2-55 
ja O336 O151 TO O-84 6-98 O94 
Sb 0-336 O151 6-70 O-85 6-66 0-96 
tia O-8238 0-139 704 2-21 6-74 2-45 
ib O825 0140 6-81 2-21 6-68 2°51 
tic 0323 0-140 6-64 2-21 6-58 2-50 
7a 0-336 0-150 6-91 0-82 6-97 0-99 
7b O337 0149 6-88 0-81 6-97 0-99 
Sa O312  O7121 7:10 4°51 705 4-98 
sb O313 O-121 7°20 4°53 6-96 5-07 
Se 0313 O7121 7-42 4°53 705 5-09 
Ya 0-300 0-104 4-03 4-47 3-78 4-83 
vb 0-301 0-104 3-96 4°25 3-67 4-58 
Ve O-301 0-104 4-00 4-24 3-69 4-88 





* Ib. UO(NO,)y per Ib. solution. 


** Runs la-2c made with aqueous feed concentration of 
0-346, remainder with concentration 0-347 Ib. UOI(NO 3) 
per lb. of solution. Inlet solvent concentration 0-049 Ib. 
UO ANO,) per Ib. of solution for all runs. 


work has shown that the distribution coeftlicient 
of uranyl nitrate between water and two typical 
extractants, diethyl ether and methyl iso-buty! 
ketone, is increased by a decrease of temperature 
[21]. However, the temperature coeflicient is 
small and the equilibrium data at 20°C, given in 
Table I [20], were used throughout the caleula- 
tions. 

The formation of a light yellow solid in the 
L).B.C. extract phase on standing overnight was 


attributed to the presence of peroxide in the 
Initially the D.B.C. was freed from 
peroxide and it must therefore have formed whilst 


solvent. 
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Table 4. Calculated data 
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Dibutyl carbitol extraction. 





Mean 
eatract 


Mean 
raffinate 
rate 
R 


me 
lb. hr 


Run 
No. 


rate 
E,, 
Ib. he 


0-064 
0-063 
0-065 
145 
O-150 
O148 
216 
220 
223 
O-s00 
O-306 
0-081 
os! 
190 
190 
OBIS 
Osl2 
0 448 
451 
0450 
O0T5 
O-O78 
O154 
0-160 
O248 
254 
O-250 
O-366 
O-362 


3-22 
3:13 
3-19 
3-22 
3-23 
3-18 
3-09 
3-18 
B15 
310 
3-12 
13-87 
13-95 
13-88 
13-93 
14°23 
13-60 
13-79 
13-87 
13-74 


0-64 
0-64 
0-64 
1-79 
1-83 
1-83 
+19 
#22 
$25 


5-26 


Ob 
10a 
lob 
lla 
lib 
lle 
I2a 
12b 


oo oF Se OF oe Ge 
zx ou 


— — 6 8 th +) =) 
_s ss ze Ss 
or Gt Se em te 


ww oo ce oe 
a a ee 


- 


- ee eee =— = Ss 


UOANOg), In 


(NTU on (ATU yg COKNOg)y Out 


O-V7 
1-01 
1-00 
0-93 
O-v2 
O-g2 
0-96 
ove 
0-93 
0-04 
0-93 
10) 
1-00 
0-08 
0-98 
0-04 
1-2 
0-98 
0-97 
o-ov 
O-o4 
o-g2 
O-O8 
ov. 
o-oo 
O-9S 
oo 
101 
11 


—-— tf tS tS 
—— 8 ts 
e ~ 
owas 


a ee en ee 
xo. 


“ww S&S & 


34 
35 
3 
‘ 
2 
2 
1 
22 


-) =) Ge Ge =2 =) =) =) =) 





* Calculated on extract phase. 


the solvent remained in the feed tank between the 
various runs. UO, 
and was probably UCO,2H,O (theoretical U,O, 
content 83-0°,), 
contamination of the sample with unchanged 


The solid contained 78-2° 
the discrepancy being due to 


This solid has been 


ether 


uranyl nitrate and solvent. 
observed in’ the 
systems [4, 8], the reaction being : 


uranyl nitrate-diethy! 


UO,(NO,), +3H,0 +0 -UO,-2H,O + 2HNO, (13) 

This reaction effectively removed the solvent 
peroxide from the extract when 
returned to the feed tank the recovered solvent 


phase and 


was free of peroxide. Several days after the 
run the unused solvent 


The slow 


completion of the last 
0.05” 
formation of hydrated uranyl peroxide is accom- 
panied by the liberation of nitric acid, which 


contained peroxide oxygen. 


o 


acts as a “ salting-out ” agent and increases the 
The amount of 
nitric acid formed was small and generally an 
insignificant change in the distribution properties 
of the recovered feed solution was observed. 
Peroxides do not form with M.L.B.K. and with 
this solvent both phases remained stable through- 
out the work. After several weeks the ketone solu- 
tion of uranyl! nitrate developed a light orange 


extraction of uranyl nitrate. 
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Table 5. Calculated data 


Methyl iso-butyl ketone eatraction. 





Mean Mean 
raffinate extract 
Run rate rate N,* 

No. R,, E,, Ib. hi (Ar), | (Amin 

lb. hr. Ib. hr. 
la 378 ORS 0-093 O-0Zs8 OAS 
Ib 3-79 OBS 0-004 0-027 O-O4 
le 3-78 OBS og 0-027 oor 
2a 11-71 O-s82 0-004 O-020 Oost 
2b 11-85 O-s4 O-101 0-020 OOS! 
2c 12-01 OSS O-O9T oz O-033 
bu 370 251 0-200 O56 OSI 
Bb sl 2-32 202 O-037 0-052 
+a 12°06 2:37 274 oz 0-044 
th 12-28 250 283 O28 OOS 
kk 12-06 2°37 276 oor O44 
via 710) Oso ool oo O46 
Sb TOO ool O-103 oo O46 
Ga 6-89 2-33 O-233 O-035 Oo 
6b O75 2.36 245 O82 O-O-4S 
tie ool 236 Ovz4 Ee 4 O-O48 
a oo oOo) O19 ono O-O4T 
Tb 6-93 0-90 0-108 0-082 0-048 
aa 7-08 475 O-382 O-038 0-054 
&b 7-08 sO OB 0-030 O55 
Sc 724 $81 O-304 0-039 0-055 
ha B01 +65 O-283 OOF 0-055 
9b 3-82 4-57 0-300 O-O41 055 
% i) 56 Osoo Oot! W055 


UOANOs)y In 


Ky Kp (HTU)oR (UTU)gp UCOYNOg)g Out 
7-12 12 5 1-51 0-99 
745 78 ts 14 o-oo 
7-27 bis 26 149 1-01 
lo 4s O56 14M 1-00 
10-78 Ho oOo ooo O-09 
o-41 28 140 1-08 1-00 
11-86 S38 256 2-26 1-00 
11-6 So 2-68 2-29 0-97 
20-20 Libchl SS 1-46 0O-99 
21-59 1b 66 1-39 101 
20-54 40 ST 1-45 0-99 
74 tou 7-71 1-56 0-99 
7 50 b7s 7:56 1-56 0-99 
15-07 10-16 375 1-88 1-01 
lies 1O-s2 tO 1-79 0-98 
1619 1-77 ss) 1-80 0-97 
77 oe 7-34 1-50 0-98 
7-22 bSI 7 86 1-35 O-97 
21-48 15-12 270 2-58 O-96 
21-41 119 2-71 2-59 0-97 
21-58 19-31 2°75 2-58 0-99 
14-75 loo 2-17 SAT 0-99 
15-3 11-6 2-00 $21 0-98 
15s 116 20r 320 0-99 





* Calculated on extract phase 


colour from oxidation of the solvent. but this 
did not occur during the experimental work 


CORRELATION OF Data 


According to the two film theory, the height of 
an overall transfer unit can be split into two 
individual film HTU’s, each of which is likely to 
be a different function of one or both phase flow 
rates. Previous experimental work in laboratory 
apparatus has shown that the height of an overall 
transfer unit in liquid-liquid extraction may be 
expressed as a power function of the ratio of the 
phase flow rates [18]. These correlations only 
apply over the narrow range of flow rates studied 
for each system. 
flow rates was possible in the dise column, it was 


involving power 


considered — that 
functions of the ratio of the tlow rates would 


expresst ms 


Since only a limited range of 


For the D.B.C. 
series the ratios of maximum to minimum flow 
rates are 45 and 8-7 and for the M.LB.K. series 
the ratios are 3-3 and 5-9 for ratlinate and extract 


correlate the experimental data. 


phases respectiy ely , 


Table 6. Disc column dimensions. 





Number of dises ; ‘ 79 
Average disc diameter 00-0484 ft. 
Average disc thickness 0-0148 ft. 


0-0833 ft. 
0-122 ft. 


Tube internal diameter 
Mean perimeter of dises 


Height of dises 3-82 ft. 
Mass transfer surface * 0-468 ft.2 





* Uncorrected for solvent film thickness or for loss of 
area at points of contact, but used throughout as an 
approximation to the true liquid interfacial area [14]. 
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-o 
MEAN RAFFINATE RATE R,, 


Variation of (1TU) jp, with raflinate rate for 
Extract rate para 
$22 Ib. hr.; 


Fig. 5. 
D.B.C. extraction of uranyl nitrate 
0-63 Ibo hr.: 2. 1-77 1b. hr. 3 


meters: 1, 2. 
4. 515 tb. her 


(a) Dibutyl Carbitol Extraction 


If the data in Table 4 are plotted in the form 
of (HTU)o» rate R 


raflinate % 
lb. hr. on log-log co-ordinates, with mean extract 


versus mean 
rate E, tb. hr. for the various series as para- 
meters, a series of parallel straight lines is obtained 
in Fig. 5. The slope of the lines is 0-76 and the 
data is correlated by an equation of the form ; 


(HTU) oy = 0 R97 BE, e. (14) 


The constants 6 and ¢ are obtained by plotting 


(HTU ox 
R,,° 76 


versus E. on log loy 


co-ordinates, which gives 


(HTU),, = 68 
By a similar procedure, 


(HTU)og = 97 (16) 


E 
where equations (15) and (16) relate to the range 


of liquor rates studied. The agreement of the 
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0°76 


(HT), = 68 ee 

















Variation of (2 7TU op and (11TU joy with mean 
extraction or uranyl nitrate, 


big. 6 


liquor rates for DBA 


experimental data of Table 3 with the above 
equations is shown in Fig. 6. 

The overall mass transfer coetlicient relation- 
ships may be obtained by a similar procedure, or 


more simply as follows : 


Since (HTU),» Mp mean) (8) 
Ky 

x R,, 

"06-122 (HTU jo 


1-20 R_°** EF O76, 


0-85 KL°* Eo 7, 


(17) 
Also Ky (18) 


where equations (17) and (18) apply to the range 
of liquor rates studied, Fig.7 shows the agreement 
of the experimental data with the above equations. 
(b) Methyl Iso-butyl Ketone Extraction 
The (HW TU), data in Table 5 are plotted versus 
™on log-log co-ordinates in Fig. 8, which gives 


the expressions, 


(HTU),,. 19) 
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R -0.36 
(HTU),» = 2-82 Fa : (20) 


r ) 
and from these equations are obtained the overall 
mass transfer coefficient relationships : 
4 . 0.36 , O64 ‘ 
K, = 421 B,,°%° BE, (21) 
4 


Kp = 2-91 R09 E66 
y m m . 


cA Fig. 9 shows the experimental data of Table 5 
t , plotted in the form indicated by equations (21) 
J Fs 
A 
¢ 


and (22). Greater scatter of the data for the 





M.I.B.K. series compared with the D.B.C. series 


arises from the errors in determining N, for the 
024 0-76 / 
K,=o-ssR,. E,. former series, 





The overall mass transfer coetlicient correlations 
obtained for the conditions studied show that both 
K,, and Ky increase with increasing flow rate 











of either phase. The extract rate influences the 
value of the transfer coeflicients to a greater 
extent than does the raflinate rate, as seen from 
Effect of mean liquor rates on Kp and Ky for the exponents in these equations. Under the same 
D.B C. extraction of uranyl nitrate. conditions of phase flow rates, A, is always 





-0°36 
(TU = 2-82 | 
Of E. 

















Fig. 8. Variation of (MTU)p, and (HTU jo, with mean 
liquor rates for M.LB.K. extraction of uranyl nitrate. 
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(r,)°"* (e m “64 


Effect of mean liquor rates on Ay and Ay» for 
M.L.B.K. extraction of uranyl nitrate 


Fig. 9 


greater than Ay, and this is because the slope 
of the equilibrium curve ts greater than unity. 
Higher overall mass transfer coellicients were 
obtained for given tlow rates with M.LB UK. com- 
pared with D.B.C. as the extractant. This is 
caused by the greater uranyl nitrate diffusivity 
in the former solvent and its lower viscosity (the 
pure solvent viscosities at 23°C are 0.57 and 2-33 
centipoises for M.1L.B.K. and D.B.C. respectively). 
Mexpocn and Prarr [9] have shown that the 
ratio of the Schmidt groups for the D.B.C. and 
M.I.B.K. phases is approximately 15. 

Figs. 6 and 7 show that (HTU),,,. increases 


rapidly as the ratio |™ is increased whereas 


(HTU),, 


increased. 


“™m 


decreases slowly as this ratio is 


(c) Individual Film Coefficients 
If it is assumed that the individual film resist- 


ances are additive and that the interfacial resist- 
ance is negligible, then 
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Rn (HTU)». 


(ATU), aE 


(HTU), + 
and 


(HTU),, HE, 


(HTU), (HT R),. 


Hence 
(ATU), — R,, 
(HTU oe HE, 


From equations (15) and (16) for the D.B.C. 
seTICs, 
(HTU),, O70R, 


; (26) 
(HTU jo E 


Therefore. the of concentrations 
used, from equations (25) and (26), the average 
In an 


over range 
slope of the equilibrium curve is 1-438, 
attempt to separate the individual film (HTU)s’, 
the data in Table 4 were plotted according to 
equation (23). Satisfactory correlation was not 
obtained and is probably due to the large curva- 
ture of the equilibrium line, preventing the use 
of a mean value of the slope “ H.” However, for 
the M.L.B.K. series, the portion of the equilibrium 
eurve is involved is almost linear and from 
equations (19) and (20), the mean value of “ H” 
is 1-45. The data in Table 5, 
shown in Fig. 10, give average values of (HTU), 
of 15ft. and (HTU), of 1-0 ft. Thus the 
individual film resistances for this system are 


the 


when plotted as 


approximately equal, which conforms with 











r 4 4 r. i 


4 ie) 





‘ 
Rm 
E.. 


Individual film 17U°s for M.LB.K. extraction 
of uranyl nitrate. 


Fig. 10. 
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results of Murpocn and Pratt [9]. Insuflicient 
information was available from this work to 
determine the magnitude of the interfacial 
resistance which is present when uranyl nitrate 
transfers from aqueous solution into organic 


solvents. 


Disc SSION 


In the runs with the lowest extract rates there 
was a tendency towards non-uniform solvent 
flow over the dise assembly, particularly if the 
column was not vertical. The solvent tended to 
low only over the centre of certain dises, but in 
general this was overcome by temporarily in- 
creasing the rate. For the runs with intermediate 
extract rates an even film of solvent was readily 
obtained. However, for those runs with the 
highest extract rates a distinct, wave-like flow 
pattern occurred over the disc assembly and was 
more pronounced as the raflinate rate increased, 
Solvent droplets occasionally broke away from 
the top 8-15 dises and these droplets were not 
collected by the arm above the top disc. With 
higher extract rates this effect became more 
pronounced and made it impossible to operate 
the column. The upper extract rates used in 
this work are therefore near the limit of satis- 
factory operation of the column and the condition 
corresponds to the transition point that occurs 
before * flooding ~ in spray or packed columns. 

The relationships for overall mass transfer 
coeflicients and heights of overall transfer units 
obtained from this study of the dise column are 
similar in form to those obtained in the extraction 
of organic solutes in either spray or packed 
columns [18]. In the dise column the dispersed 
phase wets the packing, whereas in packed 
towers the packing usually is wet by the continuous 
phase. This is an important difference in the 
operation of the two pieces of equipment and 
further work is required to determine if the 
results from dise column studies of liquid-liquid 
extraction processes can be applied to the 
prediction of operation of packed towers. 

The dise column has been shown to operate 
satisfactorily for liquid-liquid extraction, but its 
use is limited to comparatively low liquor rates 
for the phase flowing over the discs. Despite this 
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disadvantage, the use of the dise column allows 
mass transfer coefficients to be determined on an 
area basis. For solvent extraction studies, this 
column is likely to prove more useful than the 
wetted-wall column. In the latter apparatus, the 
wall liquid is not fully turbulent and it is difficult 
to obtain satisfactory mass transfer data for this 
phase. In the dise column the dispersed phase is 
probably mixed cach time it flows from one dise 
to another by the turbulence which may be 
expected to occur at the dise junctions. 
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NOTATION 
(a) Symbols 
a — surface of contact between phases, ft.? /ft*. 
A & surface of contact between phases (total area 
of dises), ft®. 
KE & extract or solvent phase flow rate, Ib. hr. 
il — slope of the equilibrium curve. 
HTU — height of individual film transfer unit, ft. 
height of overall transfer unit, ft. 

K — overall mass transfer coefficient, Ib. uranyl 
nitrate (hr.) (ft.*) (Ib. urany! nitrate Tb. solu- 
tion). 

Ny rate of uranyl nitrate transfer between phases, 

Ib. hr. 


(NTU), number of overall transfer units. 
R — raflinate or aqueous phase flow rate, Ib. hr. 
W wetted perimeter of dises, ft. 
‘ uranyl nitrate concentration in raflinate or 


aqueous phase, Ib. UO.(NO,), Ib. solution. 
(AD)om (a Leem true mean concentration driving 
force for raffinate phase. 
y = uranyl nitrate concentration in extract or sol- 
vent phase, Ib.UO.(NO ,),/ Ib. solution. 
(AWem = Ye — Dim = true mean concentration driving 
force for extract phase. 
Z = height of dises, ft. 
I’ = liquor rate, tb. hr. ft. wetted perimeter. 


bie constants. 

























(b) Subscripts 
1 = conditions at bottom of column. raflinate or aqueous phase. 


” 
- 


e 


(1) 
(2) 
(3) 


[4] 
(5) 
(6) 
(7) 
(8) 
[9] 
[10) 
[11) 
[12) 


[133] 
(14) 
[15] 
[16) 
[17] 


[18] 


{19} 
(20) 
(21) 
(22) 
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extract or solvent phase. 


= conditions at top of column. m — mean. 
at equilibrium. im logarithmic mean, 
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The paper is concerned particularly with correlation of phase compositions in 


Correlation of equilibrium compositions in quaternary two-phase liquid 
systems 





1954) 


systems in which one component is partially miscible with two others and completely miscible 
with the third. Examination of the experimental data of Prince and Hunrer [9] and CHane 


and Movurron [5] showed that : 


(1) The Brancker, Hunter, and Nasa [2] geometrical correlation did not apply generally 


to these systems, either in its entirety, or in part. 


(2) The quaternary tie lines may be related to the ternary lines by two correlations, 


requiring ternary data only. One correlation may be expressed as 


Xap | Nan 
Xpp 


K\Xpn + KoXen 


\" m (ny Ny) 


for the system A, B, C and D; B-D and C-D are partially miscible, 77 denotes the D 


poor phase, and m the ratio Neg (Xey 


X yy), where X¢jy is the wt fraction of C 


in the 7 phase. The second correlation has not been expressed algebraically. 


(3) Prediction of quaternary equilibrium required the two correlations and knowledge 


of the quaternary saturation surface ; the latter may often be approximated to, e.g. 


by CuanG and Movuvron’s interpolation. 


(4) The evidence supported CuanGc and Movuron’s plait point correlation ; but the first 
of the above correlations would make its use unnecessary, 


INTRODUCTION 

The design of solvent extraction processes and 
equipment requires amongst other data, informa- 
tion about the equilibrium states of the two 
phase liquid systems involved. While in many 
cases such systems may be treated as ternaries, 
in others, particularly those in which the solvent 
consists of several components, such as the 
“mixed” and “double ~ solvents, the actual 
system, or an adequate model, must include 
four components. If the equilibrium properties 
of such a four component system could be 
derived from those of the component ternaries, 
much experimental work could be avoided ; 
it would be possible to estimate the effect of 
introducing auxiliary solvents, and make the use 
of restricting assumptions, such as the complete 
immiscibility of two solvents, unnecessary in 
many cases, 
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The relation between quaternary and ternary 
equilibrium will therefore be examined, from an 
empirical point of view, for those systems for 
which the necessary data are available. While 
some studies of this type have been made for 
quaternary systems in which only one pair of 
components is partially miscible, less is known 
about those in which one component is partially 
miscible with two others: detailed consideration 
will therefore be given to systems of the second 
type. 

It should be noted that in a quaternary system, 
the determination of the composition of a phase 
in equilibrium with a given phase (or the deter- 
mination of the composition of two equilibrium 
phases obtained from a known heterogeneous 
mixture) involves three unknowns; and in this 
sense, a quaternary tie line has three degrees of 


freedom. Complete determination of a tie 
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line requires thus knowledge of three restraints 
to which the lines must be subject. One of these 
may be particularly associated with the saturation 
surface. and this problem will be considered 
separately from that of the other two, to which 
the earlier, and major, part of the following 


discussion will be devoted. 


PUBLISHED QUATERNARY CORRELATIONS 


On the basis of their results for the system 
acetic acid — chloroform — acetone — water at 25 C, 
Brancker, Hunrer and Nasu [2, 3] proposed 
of the 
tetrahedral model used to represent such systems 
the ternary 


triangles ; 


a correlation in terms of the geometry 


representation of 
Fig. 1). Quaternary 


(analogous to 


systems on 


Fig. 1. Single binary quaternary 


tie lines such as OP were found to lie in planes 
containing a ternary tie (MN) the 
opposite vertex (4); and they could therefore 
be defined by the intersection of such 
planes (AMN, DKL). based on the two ternary 
systems (BCD, ABC : 
acetic and acetone 

water). Further, they found the 
saturation surface to be the locus of straight 
lines joining the ternary saturation lines and lying 


line and 


two 
in their case. chloroform = 


chloroform 
quaternary 


water acid 
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in certain planes; it is not clear whether these were 
planes through the acetone — acetic acid edge 
of the tetrahedron, or perpendicular to the 
chloroform — water but the difference 
between the two may normally be neglected. 
The equilibrium properties of the quaternary 
system were then completely defined in’ terms 


edge. 


of those of the component ternaries. 

The only other detailed study of a similar 
system was undertaken by Prarr and GLover 
[8] who reported that the system. acetone — 


correlation : 


acetaldehyde acetate — water. not 
fit the but 
Braxcker [1] claimed that their calculations 


were in error, and that the system would fit 


vinyl 


above geometrical 


if these were carried out correctly. 

The correlation, as proposed, applies only to 
quaternary systems in which one pair only of 
components is partially miscible (hereafter referred 
to as “ single binary quaternaries "). In the case 
of more complex systems, e.g. those in which one 
component is partially miscible with two others 
(which will be referred to as “ double binary 
quaternaries ’), where there are three ternaries 
with two-phase areas, we would not. a priori, 
expect the apply 
with respect to all three ternaries, as any two 
could then be used to predict the third. It may, 


geometrical correlation to 


A 


<= 


~Yy 
ry 


Fig. 2. Double binary quaternary. 
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however, apply with respect to only one or two 
of the ternaries. Fig. 2 is a tetrahedral represen- 
B, and C denote 
the same components as in the two systems 


tation of such a system; A, 


discussed in detail below, and D corresponds to 
water (W) in the latter. 

A detailed study of a double binary system, 
benzene — ethyl isovalerate —ethy! alcohol — water. 
has been reported by Cuanc and Mounron [5] ; 
from their results for this, and several similar 
they that the 
correlation applied with respect to the band 


systems, deduced geometrical 
ternary (the ternary containing the three non- 
consolute components; BCD in Fig. 2). They 
claimed that the saturation 
surface could be obtained by a suitable straight 


also quaternary 


line interpolation between the ternary saturation 
lines, and that the quaternary plait points could 
be found from the ternary plait points. However, 
even if the three correlations involved were 
Table 1. 


Benzene 
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generally valid for these systems, they would 
not be suflicient to predict quaternary equili- 
brium from ternary data only. 

There are data for three other double binary 
isobutene 
naphtha at 20°F [11] and methanol 


carbon disulphide 


systems: in the case of two, butadiene 

furfural 
oleic acid — olive oil at 25°C 
[7], at least one “ component” is a mixture of 
the third, 


cyclohexane at 25°C [4], the 


several compounds ; for aniline 


benzene — cetane 
published tie lines have been calculated making 
unnecessary assumptions, but it was found pos- 
sible to recalculate these on a more valid basis. 
Analysis of the data for these systems suggested 
that the geometrical model did not apply in 
respect of the open ternaries (ABD, ACD, in 
Fig. 2), but that it might apply to the band ter- 
nary. The specific claim that the model in fact 
applied to the band ternary was made by Moreno 


~ 


and Crespt Gonzaes [7] for their system. 


Equilibrium compositions for the ternary systems. 
Acetic Acid 


Water Ternary 





Phase 


B 


07-22 
97-20 
O4-84 
94-20 
74°77 


72-67 


W Phase 


B é Ww 


O57 
OST 
0-938 
104 
10-00) 
11-24 


73-87 
73-81 
62-58 
60-65 
24-65 


23-03 


65°35 


65-73 





Carbon Tetrachloride 


Acetic 


Acid — Water Ternary 





Phase 


99-62 
99-02 
97-93 
97-02 
05-44 
93-14 
89-90 
87-21 
83-32 


O36 
0-96 
2-04 
2-04 
449 
6-76 
9-92 
12°54 


O-o2 
0o-02 
0-03 
0-04 
0-07 
0-10 
O-18 
0-25 
0-42 


W Phase 


15-12 
26-97 
40-91 
48-02 
56-22 
62-53 
66°39 
67°27 


66-70 
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Table 2. Equilibrium compositions for the quaternary system. 





Phase W Phase 


87 tS? é 38-53 60-45 
21 24:77 2: : 66-05 20-91 
26 2:50 2 26-15 73-29 
23:71 HR ‘ 38-84 60-12 
7s 21-08 . 65-47 24°51 


O-O5 1-80 ‘2 26-09 73-37 
48-20 3:79 3 38-08 60-00 


te 
-_ 


te 
> 
2 


45-95 8:57 55-54 41-18 
42-25 17-71 } 2: 65-61 24-32 
39-52 23-00 far 66-05 18-08 
49-08 48-95 1 2: 26-37 73-08 
47-96 47-99 34 3f 30-24 59-71 
39-74 41-47 Is-14 . 65-67 24-00 
36-35 38-27 24-238 ¢ 65-85 18-12 


fo to to te te ot 
2 aude & 


= 


73-61 1-78 $ 30-28 69-09 
72-70 301 30-59 59-38 
TO-02 6-2 2 : 56-46 40-08 
64-64 14-85 66-06 23-54 
“< o1eol 19-71 ‘ HH24 17-86 
24-67 73-79 149 3: 27-00 72-46 


RSs 
tS to to to 
i 
c— © te oO 
ae ae es | 


ts 
+ 
= 
x 
x 


24-08 72-80 B05 39-08 58-95 
20-19 63-63 15-73 3-7: 77: 66-29 22-21 


19-00 59-63 20-60 : 66-01 17-26 


=* 
~ 
~ 
= 





To test the above correlations, and explore ACETIC ACID 

the extent to which the ternary data may be 
used to predict quaternary equilibrium, tie line 
data were obtained for the system benzene 

carbon tetrachloride — acetic acid — water at 25°C 
by methods to be described elsewhere [9]; 
they are shown in Tables 1 and 2. The relation 
of the quaternary to the ternary tie lines will 


now be discussed. 


THe GromerricaL CORRELATION 


The application of the correlation to this system 


was examined by calculating the lines of inter- 
section of planes through a quaternary tie line 
and a tetrahedral apex with the opposite ternary 
plane, for the three ternary planes containing — seNzene RS I 
solubility gaps; if the correlation applies, these TERMARY THE LINES 
lines should be tie lines of the corresponding 

ternary systems. Such lines of intersection Fig. 3 
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(which will be referred to as “ projected tie 
lines’), found for a representative group of 
quaternary tie lines, are shown in Table 3; 
for convenience, the lines are characterised by 
their intersections with the sides of the conven- 
tional triangular diagrams (edges of the 
tetrahedron), given as weight °, of the appro- 
priate component. Figs. 3, 4 and 5 show the lines 


ACETIC ACID 








CARBON TETRACHLORIDE WATER 


Fig. 4 
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plotted on the triangular diagrams ; the quater- 
nary lines being identified by the numbers 
given to them in Tables 2 and 3. 

In the two open systems (benzene — acetic 
acid — water. Fig. 3; carbon tetrachloride — 
acetic acid — water, Fig. 4) the projected lines 
depart progressively from the ternary tie lines as 
the quaternary compositions depart from the 


WATER 












18 1920 17 26 27 23 28 29 30323837 


BENZENE CARBON TETRACHLORIDE 
Fig. 5 


Table 8. Projection lines for the system Benzene (B)— Carbon Tetrachloride (C) — Acetic Acid (A) —- Water (W) : intercepts 
with binary lines as weight ©. 





Line Plane ABW Plane ACW Plane BCW 
Origin A Origin A Origin B 
W B W Cc Ww Cc 

7 90-04 26-01 
18 73-71 96-87 73-71 91-15 99-89 24-88 
19 60-79 03-78 60-78 83°31 YS -S4 24-94 
20 28-05 75-63 27-82 51-34 91-84 25-38 
2: 97-68 50-39 
26 73-49 96-32 73:49 96-31 99-93 49-93 
27 60-38 92-36 60-36 92-54 99-55 50-02 
28 27-53 70-75 27-28 71-65 92-65 51-10 
29 22-04 64-09 22-60 65°35 89-20 51-38 
30 99-92 74-98 
32 98-57 75-34 
35 72-86 94-37 72-86 98-05 99-90 74-94 
4 59.62 88-01 59-61 96-04 99-74 75-14 
37 25-88 58°30 25-60 81-52 94-49 75-94 





22-09 21-60 





7701 90-28 75-92 
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ternary, and some of the projected lines even 
intersect each other, In the band = system 
(benzene — carbon tetrachloride water), the 
mutual solubilities in the two binaries are very 
low (the highest is that of benzene in water, 
O18 g per 100 g), so that the water phase consists 
of at least 99-8°,, water, and tie lines may be 
regarded as passing through the water apex. 
If the geometrical correlation then applies to this 
ternary, quaternary tie lines must lie in planes 
through the water and acetic acid apeces, t.c. planes 
of constant ratio benzene : carbon tetrachloride, 
and this ratio must therefore be the same both for 
phases. Table 2. however, shows marked and 
systematic differences in this ratio; these devia- 
tions between ternary and projected tie lines 
may also be seen in Fig. 5. It may be noted 
here that as benzene and carbon tetrachloride 
have very different densities (0-87, 159g mil.), 
and as essentially the density was used to deter- 
mine the ratio of these components in the water- 
poor phase, these figures were determined com- 
paratively accurately. 

For all ‘three ternaries the deviations are thus 
outside the expected limits of experimental 
error, and from their systematic nature, it would 
appear that the geometrical model does not 
apply to the system, and that here planes through 
ternary tie lines and the opposite vertex do not 
contain any quaternary tie lines. 


Provosep Corretarion | 


To obtain a more suitable correlation for double 
binary quaternaries, it was decided to modify 
the logarithmic correlation of Haxp [6], this 
being the most generally valid of the ternary 
tie line correlations. For ternaries, Hanp 


proposed the relationship 


Nas Ney — RUN yr Nyy)" (1) 
where 

weight fraction of (consolute) com- 
ponent R in the S rich phase ; 


X ps 


K, » = constants for a particular system. 
In a ternary system, a straight line should then 
be obtained if the above fractions are plotted 
on a logarithmic scale. 
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Denoting, in the system studied, the combined 
benzene and carbon tetrachloride in any phase 
as the Hf 
of this H “ component” as weight 
tetrachloride, by M, it was found, on plotting 
log (XN gy Xipy) against log (XY ,,, X,,,,). that 
the two lines representing the two open ternaries 
were parallel, and that points representing 
quarternary lines of constant HW) phase MV lay 


“component,” and the composition 


» carbon 


on straight lines, parallel to the ternary lines. 
The relationship then held for the system : 


XN wi Nan A’ (X in“ Nw " (2) 


mn being a constant, and A’ dependent only on the 
composition of the HH component of the H 
rich phase. The relation between A’ and the 
composition was found from the spacing of the 
constant M lines; and hence it was found, 


that for the system 


XN 4H Nyon A(X au Niu y" (3) 
where 

Nyon Xun k + Xun 
K, », k ~ constants for the system: 0-1000, 


1-214, 3-10, resp. 

Fig. 6 shows a plot of equation (3) and of the 
results for the system. Due to the coordinates 
used, it represents an ordinary Hand plot for 
the benzene ternary; the Hand line for the 
carbon tetrachloride ternary has been superim- 
posed on the benzene line by use of the fictitious 
H concentration Xj 4,, The extent to which 
equation (3) applies to the quaternary tie lines 
may now be gauged from the deviations of the 
corresponding points from the line ; such devia- 
tions are quite small. Apparently systematic 
deviations occur at the highest acid concen- 
trations ; the effect is probably that often found 
on ternary Hand plots, as the plait point is 
being approached, and is due partly to the greater 
experimental error in the tie line determinations, 
and partly due to systematic errors in the satura- 
tion curve. The experimental uncertainty in 
this region is comparatively large: as here X yyy 
is relatively small, such uncertainties have a 
relatively large effect on the correlation. 

It is clear from the form of equation (3) that 
the three constants involved may be (and in the 
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present case were) determined from ternary tie 
line data only; here, A and » are constants in 
the corresponding Hand equation for the benzene 
ternary, while & is obtained by dividing A’ for 
the benzene ternary by A for the carbon tetrach- 
loride ternary. 

The system used to derive equation (3), and 
hence the equation, represents a special case, 
in that generally, the two open ternary systems 
We would 


the general 


would not have parallel Hand lines. 
expect, however, 
similar equation would apply if the constant ” 


that in case a 
were replaced by a suitable function of the com- 
position, say M. This supposition may be inves- 
tigated using the data of Cuanc and Mouvron [5] 
ethyl isovalerate — ethyl 
was 


on the system benzene 


alcohol — water. As for this system, M 


reported to have the same value for both phases 
in equilibrium there is no ambiguity as to the 


Correlation I applied to the system benzene 
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carbon tetrachloride — acetic acid — water. 

appropriate value to use; and as the quaternary 
tie lines had been determined as series of con- 
stant M, the hypothesis could be tested by making 
ethyl 
The series 


an ordinary Hand plot, treating benzene 
isovalerate as one component (/7), 
of constant MW were found to give straight lines 
(within experimental error), with slopes changing 
regularly on traversing the quaternary surface, 
thus showing that the line slope was a function 
of M. Further, it appeared that the change 
in slope was proportional to M, so that the 
correlation equation could now be written : 


XxX AW 
XxX wi 
. , Xan - ) nim (ni ma) (4) 
A ix BH ~~ K,X CH. 
where 
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ethyl — aleohol, 


ethyl 


components 
benzene, isovalerate, 


water. 
water poor phase. 


M 100 wt 
benzene). 


ester (wt ester 


constants for the system : 
0.960, 0-660, 0-107, 0-731 resp. 


n,. Ny Ky. Ke 


These constants may all be derived from ternary 
data only : 


the 


ester, 


corresponding 


Hand 


slopes of 
benzene, ternary 


lines. 


1k," K of the 


Hand line. 


benzene ternary 


Fig. 7 shows a plot of equation (4) and of the 
As for Fig. 6, 
dinates have been so chosen that we 


results for the system. the coor- 


have the 





=, log Xan/4ww 
° 


nM 
a, mA 











is ro 


log mn —- 
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ethyl alcohol 


Fig. 7. 


ethyl isovalerate water. 


normal Hand plot for the benzene ternary, and 
the Hand line for the ester ternary has been 
superimposed on the benzene line (as before, 
Xn + KReXey/K,). Ut appears that 


x H*H 
the deviations of the quaternary points are of the 


same order as deviations of the ternary points, 
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and hence equation (4) holds for the system, 
within experimental error. 

Although the correlation the 
required condition that all parameters may be 


above meets 
caleulated from ternary data only, it is not 
the 
equilibrium compositions, even if used in conjunc- 


suliicient for calculation of quaternary 


tion with knowledge of the saturation surface, 


as we have only two restraints for three variables. 
Prediction therefore requires a second correlation. 


Proposep CorRELATION II 


If the data for the benzene — carbon 


tetrachloride 


system 
acetic acid — water are accepted 
as indicating that the H component composition 
differs in the two equilibrium phases (as against 
Cuanc and Mou .ron’s hypothesis that it does 
not) a correlation is required in terms of this 
composition. The ternary analogy would be 
correlation of tie lines in a band system; the 
most convenient correlation is one proposed by 
Varreresstan and Fenske [12]. The quantity 
to their 8, 
X pw Ncu/A pn Xcw- denoted by @. 


@ was calculated for the quaternary tie lines ; 


corresponding in the quaternary 
will be 


and for the higher acid concentrations, found 
to give a smooth plot against H7 phase acetic acid 
concentration. This plot was related to non- 
quaternary data as follows, 

Owing to the low solubilities in the benzene 
carbon tetrachloride — water system, no experi- 
mental ternary data are available : but Q (or 8) 
may be predicted for the system from binary 
solubilities, if we assume that the benzene and 

activity 
water 


tetrachloride coellicients 
the 


they are independent of the benzene : 


carbon are 


dependent on concentration only, 
that is, 
carbon tetrachloride ratio. As these components 
almost 


form ideal mixtures, this assumption 


should not lead to great error. Now, owing to 
the low solubilities in the system, the water 
either 
the 
tetrachloride coeflicients are also constant. As 


concentration in phase is effectively 


constant. and hence benzene and carbon 


the activity of any species must be the same in 
the two equilibrium phases, we have : 


“aw Yaw ~ “py Yaw 
“nw “pn Yau Yaw Ky 
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and K, may be found from binary data. Then 
R= Xaw Xen Xan Xen 
“pwton “pn cow 
K, K, 0.000414 16 
C-O00090 


where the binary data were obtained from 
SEIDELL [10}. 

The same argument should hold at very low 
concentrations of acetic acid, and two sets of 
values of Q were calculated, using experimental 
isotherm data. The logarithms of these calculated 
values were then plotted against HM phase acetic 
As @ must be unity on the 


plait line, we can locate that end, after calculation 


acid concentration. 


of the plait points as described below, and hence 
draw smooth curves through the calculated points ; 
these curves, together with the experimental 


values of Q, are shown in Fig. 8. Experimental 
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Vig. 8. Correlation Il applied to the system benzene 


carbon tetrachloride — acetic acid — water. 

points for low acetic acid concentrations have 
not been plotted; reference to Table 2 shows 
that the benzene and carbon tetrachloride con- 
centrations in the aqueous phase are so low that 
little significance can be attached to the measured 
value of Q. 
the experimental Q 
and reasonable agreement can be seen between 


At the higher acid concentrations, 


becomes more reliable, 


the experimental and predicted values. 


Piarr Poinr CorreLATION 


Correlation I and the experimental isotherm 


are sullicient to determine the quaternary plait 


Is: 
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> 
the two ternary, and three quaternary 
plait points found by this method for the system 


line : 


studied, are shown in Table 4. 








Table 4. Plait points in the quaternary system. 
M Hi B Cc A Ww 
0 43-80 43-80 49-49 6-71 
25 16-34 34-75 11-59 47-66 6-00 
ww 49-42 24-71 24-71 45-38 5°20 
75 53-08 13-27 39-81 42-65 4:27 
100 58-44 58-44 38°35 3-21 
Cuanc and Movunron found that. for the 


systems they investigated, the quaternary plait 
line lay in one plane, passing through the consolute 
component vertex (acetic acid in the present 
Hence the locus of points of inter- 
the 
and the plait points with the (opposite) base 


system). 


section of lines through consolute apex 


should be a straight line. This was tested for 
the present system, by Cuanc and Mouvrton’s 
method of calculating Wo(B C W) for a 
quaternary plait point of a given value of M, 
using the known positions of the ternary plait 
the straight 
line referred to above); and comparing § this 


points (and hence establishing 


figure with the experimental ratio, as shown 

in Table 5. The experimental figures are all 

Table 5. Cuanc and Movuvron plait point correlation : 

calculated and experimental values of W(B ~~ C + W) 
for the quaternary plait points of Table 4. 





M Calculated Experimental 


25 0-1140 0-1146 
50 00-0944 0-0951 
75 00-0736 00-0745 





higher than those predicted, but the differences 
are small, and could be taken as due to experi- 
for the saturation surface 


mental error, which, 


in this region, is comparatively large. 


SATURATION SURFACE CORRELATION 
Cuanc and Mouiron showed, for their system, 


that the quaternary saturation surface could be 
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approximated by straight line interpolation 
between ternary points having the same ratio 
of water to H * component.” agreement between 
the experimental and interpolated values being 
within 5°,. 

In the present system, the two ternarics have 
that it is diflieutt 
to test any particular method of interpolation : 
but the the 


(see [9]) suggests that straight line interpolation 


very similar solubilities, so 


appearance of solubility curves 


would result in systematic deviations Some 
calculations were made to determine the order 
of these value of M 
and of WH (which between them locate a point 
the 


caleulated by 


deviations : for a given 


on the quaternary surface), acetic acid 


concentration was interpolation 


and compared with the experimental value 


these figures for three points are shown in Table 6, 


Table 6. CHanc and Movwron isotherm interpolation. 





Wo A cate. 1 exp 
16-78 
HORS 


GO155 


OM677 
O-4284 
O-42584 


16-28 
ooo 


H2-05 





acid concentration. while 


The 


small, 


differences im 


are significant. and systematic. These 
differences may be related to the experimental 
error by calculating the corresponding change in 
the component directly determined, in this case, 
water; as shown in the last column of the table. 


The differences are somewhat greater than the 


expected errors over those regions of the satura- 


tion surface. 


Disc UL SSTON 


If the saturation surface is known, prediction of 


equilibrium in quaternary systems requires 


knowledge of two restraints, and two correla- 
tions have therefore been proposed. applicable 
binary systems, 


to ~~ double quaternary 


Correlation I be expressed algebraically, 


and the four constants involved may be obtained 


may 


from ternary data only four accurately 


determined ternary tie lines are sullicient for 


this purpose. From the general validity of the 


Ist 
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corresponding ternary equation, and the good 
agreement between experimental and calculated 
results in the two systems for which suflicient 
the should be 
applicable to double binary quater- 


data are available. correlation 
generally 
naries; although there may be some doubt 
as to the appropriate value of m to use in equation 
(4); as this equation has been based on a system 
for which m had the same value in both phases, 
In the general case, we would expect that the 
that of the H 
phase, as that was the relevant value in equation 
(33). 


expect. m differed somewhat in the two phases 


correct value to use would be 


It may also be noted that if, as we would 
in Cuanc and Mouvron’s system, their reported 


been close to that of the H 
phase, as that phase contains the major part of 


value would have 


the relevant components : and as the reported 
value fitted equation (4), the true H phase value 
should also fit. 

Correlation IL ts not as satisfactory, in that it 
is only represented by a series of curves, and 
that for the svstem for which it was established, 
these curves could not be found very reliably. 
This was largely due to the very low solubilities 
in the band ternary, as a result of which no 
values of Q could be 
obtained near that end; also, the consequent 

addition of 
acid) 


rehable experimental 


rapid increase im solubility on 
caused 


The 


results of Cuanc and Moutron show that. to a 


consolute component — (acetic 


marked curvature of the correlation lines. 


fair degree of approximation, this correlation may 
often be replaced by the assumption Q = 1. 
Kven if agreement between correlation II 
and the experimental results may not be very 
satisfactory quantitatively, the qualitative agree- 
ment supports the conclusion from the experi- 
mental figures, that @ is not unity for the system 
involved, and hence that the Bkrancker, HUNTER, 
and Nasu geometrical correlation applies to these 
systems neither in respect of the band ternary, 
nor of the open ternaries, and therefore that the 
the 


suggested by by Cranc and Mouvron, is only 


restricted application of correlation. as 


an approximation, 


CHanc and Movunron’s correlation of the 


saturation surface may also be an approximation, 
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although here the evidence is not as clear; 
the use of this correlation should lead to com- 
paratively little error, particularly if the two 
open ternaries have similar solubilities. 

The results for the 
system fitting CHANG 


experimental present 


can be taken as and 
plait but it 


should be noticed, that in our case, the quaternary 


Mouron’s point correlation : 
saturation surface in the plait point region is 
almost parallel to the correlation plane, so that 
a decision as to whether the correlation was 
applicable would require much more accurate 
data. Further, if there were any real deviation 
from of the 


consolute would be 


the correlation. the concentration 


component in’ particular 
subject to comparatively great error. The use 
of correlation I, which requires no expermental 
plait point data, and for which other evidence is 
available, seems preferable in locating plait points. 
Quaternary tie lines may then be predicted 
in these systems by using correlation I; correla- 
tion Il or the approximation Q — 1; and the 
When the 


latter may be estimated by interpolation between 


experimental saturation surface. 


the ternary saturation curves, no quaternary 


data are required for complete prediction of 
equilibrium. It should be added, 
that 


may be derived entirely from data on the two 


quaternary 


however, quaternary correlations which 
cannot allow 
the 


not common to those ternaries, The pair involved 


open ternaries for the specitic 


interaction between pair of components 


is the pair of components partially miscible 
with a third component, and in’ the systems 
discussed above, we would expect this approxi- 
mation to lead to negligible error, and for systems 
of this type generally, the errors should be small 
compared to other uncertainties. Even correla- 
tions based on the three ternaries must always 
represent approximations, although the errors 
would be of a lower order, and therefore even 
less likely to become significant. 


CONCLUSION 


The evidence discussed above suggests that in a 
quaternary system of (liquid) components 
A, B.C and D, where B-D and C-D are partially 
miscible, the relation 
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mim (mi — M2) 


Xap 
Xpp KiXcu 
holds: H denotes the D poor phase, and m 
the ratio X,4, (Ney + X yy)- 


Xuan 
| KX py 


A second correlation, involving particularly 
the concentrations of B and C in the two phases, 
is also suggested. 

The only quaternary data required when using 
these correlations in the prediction of quaternary 
equilibrium are then saturation data, which are 
comparatively easily obtained experimentally, or 
may often be found with little error by inter- 
polation between the ternary saturation curves. 

The Brancker, Hunter and Nasu geometrical 
correlation does not appear to apply generally 
to these systems, cither in its entirety, or in part. 

While the Cuaxc and Movunron plait point 
correlation may apply, plait points should pre- 
the first correlation 


ferably be estimated by 


proposed, 
The 


largely carricd out in the Chemical Engineering 
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NOTATION 
Component 1; acetic acid or ethyl aleohol, as 
appropriate ; wt °,, of A. 
°, of B. 
tetrachloride 
wt °, of C. 


Component B ; benzene ; wt 


Component C: carbon or ethyl 


isovalerate, as appropriate ; 
Component D. 

il * component.” B ‘; Band C — rich phase, 
i.e. D poor phase. 

Fictitious 77) component concentration required 
for correlation 1. 
Constant, particularly in 


Constant. 


HAND'’s equation. 


Composition of 17 component, C H/, as wt wt °,. 
M 100, 
Constant in HANb’s equation. 
New Xen Xpu Xew- 
Ww Water: wt °, water. 
Xps — Wt fraction of component K in the S rich phase. 
'rs — Mol fraction of component #& in the S rich phase. 
yrs — Activity coefficient of component RK in the 8 
rich phase. 
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ANNOUNCEMENT 





Oxidation Conference —1954 


A joint conference of the institution of Chemical Engineers and the 
Chemical Engineering Group of the Society of Chemical Industry in 
London with the Royal Institution of Engineers (Chemical Engineering 
Group) of the Netherlands and the Royal Netherlands Chemical 
Society (Section for Chemical Technology) on Oxidation Processes 
was held at Amsterdam on May 6th and 7th, 1954. 


The publishers (Pergamon Press) have in preparation the proceedings 
of this conference which they intend to publish as a supplement to 
Volume 3 of CHEMICAL ENGINEERING SCIENCE by about 
October or November of this year. Details of price, size and a firm 
date of availability, will be announced later. 


The following papers read at the conference will be published in the 


supplement together with the relevant discussion. 


Prof. Dr. H. I. Waterman 


Introduction to the Conference 


Dr. P. Mars and Prof. Dr. D. W. van Kreve_en 


Oxidations carried out by means of vanadium oxide catalysts. 


Herrn. Dr. G. Wrerzen 


Herstellung synthetischer Fettsduren durch Oxydation von 
paraflinischen Kohlenwasserstoffen mit  molekularem  Sauerstoff. 


Prof. M. B. Doxnatp and Mr. J. C. Grover 


The air oxidation of benzene to phenol. 


Dr. F. J. Dent 


The use of oxygen in gasification 


Mr. W. J. B. Cuarer 


Oxygen processes in the Steel Industry 


Dr. A. M. CLARK 


The provision of oxygen for industrial processes 
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Announcement 


Mr. R. M. Coie and Mr. A. W. Farrearrn 


The development of a process for the manufacture of p-tertiary- 
butyvlbenzoic acid, 


J. P. Forruis and Prof. Dr. H. lL. Warerman 


Preparation for cumene hydroperoxide 


Prof. E. Barruotome 


The BASF-process for production of acetylene by partial oxidation of 


gaseous hydrocarbons. 


Mr. A. F. Linpsay 


Nitric acid oxidation design in the manufacture of adipic acid from 
cyclohexanol and cyclohexanone 


Dr. G. H. Twice 


The mechanism of liquid-phase oxidation 


Prof. F. H. Garner 


Concluding Remarks 





